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NAMER - A FORTRAN IV PROGRAM FOR USE IN OPTIMIZING DESIGNS OF 
TWO-LEVEL FACTORIAL EXPERIMENTS GIVEN 
PARTIAL PRIOR INFORMATION 
by Steven M. Sidik 
Lewis Research Center 

SUMMARY 

NAMER can be used to find the Bayes procedure for designing two- level fractional 
factorial experiments given partial prior information. The required prior information 
is 

(1) A statement for each parameter giving a prior probability that it is not zero 

(2) A statement of the probability of stopping at each contemplated stopping point 

(3) A statement of the value to the experimenter of an unbiased estimate for each 

parameter 

The steps of the design and performance of the experiment may be represented as a 
finite discrete game between the experimenter and nature. The decision space E for 
the experimenter consists of the choice of initial defining parameter group, the choice of 
the sequence or sequences of subgroups that define the telescoping, the choice of 
physical- design variable matching, and the choice of parameter- estimator matching. 

The decision space N for nature consists of the choice of which of the parameters are 
zero and the choice of the stopping point of the experiment. The Bayes procedure maxi- 
mizes the expected utility over all possible distinct choices of parameter- estimator 
matchings, physical- design variable matchings, and defining parameter groups for an 
assumed strategy for nature. 

This report presents an algorithm and a computer program entitled NAMER which 
computes the expected utility of all possible physical- design variable matchings and 
parameter- estimator matchings for a specified choice of defining parameter groups. 

The matchings which maximize the expected utilities are saved and printed out. The 
computational procedure utilizes the group properties of the parameters and the stand- 
ard ordering. Complete program documentation is presented including sample input and 
output and a sample problem illustrating the usage (appendix A), program listings (ap- 
pendix B), a program symbol table (appendix C), and a general flow diagram of the com- 
puter program (appendix D). 


INTRODUCTION 


The two- level fractional factorial designs represent a class of designs of experi- 
ments which yield estimates of first- degree effects and interactions for a small amount 
of experimentation. The main disadvantage of this class of designs is that the estimates 
(using linear least- squares estimators) are always estimates of aliased combinations of 
parameters. To make conclusions about single parameters it is necessary to have some 
information about the parameters from a source other than the experiment. If such in- 
formation is available before the experiment is performed, it may be incorporated into 
the design of the experiment. 

There are many situations in practice in which an experimenter may have varying 
amounts of information concerning the variables he wishes to investigate. Sidik and 
Holms (ref. 1) have developed some optimal design procedures when the prior informa- 
tion is 

(1) A statement for each parameter giving a prior probability that it is not zero 

(2) A statement of the probability of stopping at each contemplated stopping point 

(3) A statement of the value to the experimenter of an unbiased estimate for each 

parameter 

The steps of the design and performance of the experiment may be represented as 
a finite discrete game between the experimenter and nature. The decision space E for 
the experimenter consists of the choice of initial defining parameter group, the choice 
of the sequence or sequences of subgroups that define the telescoping, the choice of 
physical- design variable matching, and the choice of parameter- estimator matching. 

The decision space N for nature consists of the choice of which parameters are zero 
and the choice of the stopping point of the experiment. 

The Bayes procedure maximizes the expected utility over all possible distinct 
choices of parameter- estimator matchings, physical- design variable matchings, and 
defining parameter groups for an assumed strategy for nature. 

This report presents an algorithm and a computer program entitled NAMER which 
computes the expected utility of all possible physical- design variable matchings and 
parameter- estimator matchings for a specified choice of defining parameter group or 
groups. The computational procedure utilizes the group properties of the parameters 
and the standard ordering of the parameters. 

The program can handle experiment designs for as many as nine factors and 32 
stopping points. The relation among the stopping points is arbitrary so that, by proper 
input of data, multiply telescoping designs may be considered or as many as 32 single- 
stage designs may be analyzed simultaneously. The program output gives the physical- 
design variable matchings and the parameter- estimator matchings which are the Bayes 
decisions. Also those matchings which maximize the expected utility at each individual 
stopping point are printed out so that a security strategy may be specified. 
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If an experimental program has already begun so that a physical- design variable 
matching is specified, NAMER may still be used to change the choices of telescoping 
options based upon revised prior probabilities of stopping at each stopping point not yet 
reached. 

The algorithm and program are fully described. Listings, sample input and output, 
and a sample problem illustrating the usage are given. 


SYMBOLS 


B 

B(h) 

h 

i®B(h) 

n 

P[A] 

Pr(A) 

Pb 

Pi 

Psh 

U 

U(h) 

U(i,k) 


u.(h) 

X A’ X B’ 

X l’ X 2’ 

® 

Y 

^I’^A’^B’ ‘ ' * 


full parameter group 

XI- 

subgroup of B used at the h Ln stopping point of experiment 

denotes stopping point of experiment 

set of standard- order subscripts of elements of /3j(x)B(h) 

number of factors (independent variables) 

permutation of ordering A 

probability of event A 

prior probability of a block effect not being zero 
probability that 0^ is not equal to zero, Pr(/3^ * 0) 

ii- 

probability experiment will stop exactly at hr n stopping point 

maximized expected utility over stopping points for a given defining 
parameter group and matching of variables 

1U 

maximized expected utility of hi stage for a given defining parameter 
group and matching of variables 

expected utility gained by assigning estimator for alias set /3j(x)B(h) 
to 0^ kei®B(h) 

xu 

utility assigned to an unbiased estimate of 0^ at h n stopping point 

independent variables (design) 

independent variables (physical) 

group operation 

dependent (response) variable 

parameters in a model equation relating design variables to dependent 
variable 
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/3j®B(h) 

A()’ * “ 


5 

e 


coset (alias set) obtained by multiplying all elements of B(h) by j3j 

parameters in a model equation relating physical variables to dependent 
variable 

random variable with mean zero and finite variance 
element of 


REVIEW OF BAYES PROCEDURE AND STATEMENT OF COMPUTING PROBLEM 


In a full factorial experiment with n independent variables X^, each re- 

stricted to assuming only two values, there are 2 n possible distinct combinations of 
values. It is common practice to say the independent variables can assume either a 
"high” level or a "low" level. Each of the 2 n distinct combinations of levels is called 
a treatment combination. From such an experiment it is possible to estimate the /3*s 
in an equation of the form 

Y = &1 + ^A X A + /3 B X B + i3 BA X B X A + /3 C X C + /3 CA X C X A + /3 CB X C X B 


+ A 


CBA X C X B X A + 


+ P 


CBA’ 


* X C X B X A + 5 


(1) 


where 6 is a random variable with mean zero and finite variance. (Note that the order- 
ing of the subscripts is the reverse of that normally used. The reason for this will be 
explained shortly. ) 

A regular fractional replicate of the full factorial design does not allow separate 
estimation of all the /3's. Certain linear combinations of them can be estimated, how- 
ever. The particular set of linear combinations which can be estimated depends upon 
the treatment combinations composing the fractional replicate or, equivalently, upon the 
choice of the design of the experiment. For example, a one- half replicate experiment 
on four independent variables would provide eight estimators which might be estimators 
of (depending upon the particular fraction): 


(Ai + ^dcba^ + %>ba) 

(/3 a + Pvqb) ^CA + %b) 

(/3 b + ^DCA^ ^CB + ^DA^ 


(2) 


(A B a + Adc^ 


(Acba + Ad^ 



From such estimators, nothing can be inferred about any single parameter without mak- 
ing some assumptions about the other parameter in the alias set. 

The set of all 2 n contrasts which provide estimators of the parameters in a full fac- 
torial form a group under the appropriate operation. There is a one-to-one mapping 
from the group of contrasts onto the group B of parameters. Since the point of view of 
this report is based upon knowledge about parameters, it is more convenient for the 
development to be in terms of the parameter group. The operation defining a group with 
respect to the parameters is analogous to that used in the group of contrasts. 

With every regular fractional replicate there is associated a defining parameter 
group (d. p. g. ) which can be used to determine the aliased sets of parameters that can 
be estimated. Conversely, given a d.p.g. , there is a regular fractional replicate asso- 
ciated with it. 

Holms (ref. 2) and Holms and Sidik (ref. 3) present a technique called telescoping 
sequences of blocks. This allows an experimenter to perform a factorial experiment in 
stages, where the starting stage is a small fractional replicate and the final stage is 
some larger fraction. Each succeeding stage adds treatment combinations to those run 
in the preceding stages. In order to retain the orthogonality and the orthogonal blocking, 
each stage must be a power of two times the size of the preceding stage and all the treat- 
ments run must form a regular fractional replicate. In what follows, we will restrict 
ourselves to single telescoping and consider the n n stopping point to be the n n stage. 

In the case of multiple telescoping, this would not be true, in general, for there could 
be many stopping points in a stage and the relations between groups are more complex. 
This is not essential to the discussion, however; and we consider single telescoping only 
to keep the notation simple. 

xi_ 

Denote the d. p. g. at the starting stage as B(l) and the d. p. g. at the h Ln stage as 
B(h). If the d.p.g. 's are such that B(h + 1) is a subgroup of B(h), the sequence of regu- 
lar fractional replicates corresponding to them will form a telescoping sequence of 

XT— 

blocks under the rules established in Holms and Sidik (ref. 3). At the n n stage, the 
treatment combinations run should form the fractional replicate defined by B(h). The 
fractional replicate at the h + 1 stage can be achieved by adding to the treatment com- 
binations defined by B(h) those treatments in the replicate defined by B(h + 1) but not yet 
performed. 

As the experiment progresses through the stages, the number of alias sets in- 
creases, while the number of /3’s in each alias set decreases. If the final stage is the 
full factorial, each parameter is separately estimable except that certain of the param- 
eters are confounded with blocks. Whether block effects physically exist is a question 
the experimenter must answer. 

It will be convenient at this point to introduce an alternate notation for equation (1). 
Let the n independent variables be denoted as X^, • • • , X n> Number the 2 n jS’s of 
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and consider the following equation which is similar to 


equation (1) from 
equation (1): 


/3 0 t0 P n 

U 2 n - 1 


Y = /3 q + /3 lXl + P 2 X 2 + + ^Xg + • • • + 0 


„ XX - 
2^1 j n n- X 


x i + 


( 3 ) 


Equation (1) and equation (3) are both written in what is called the standard order. If 
the subscripts of the 0's are rewritten as n- digit binary numbers, it becomes quite 
obvious how the terms and coefficients of equation (3) are related. For example, let 
n = 4 and consider the following equation, where the subscripts on the (i's are written 
as binary numbers: 


Y = + Vi + ^io X 2 + % X 2 X 1 + ^100 X 3 + ^101 X 3 X 1 + 

+ ^111 X 3 X 2 X 1 + ^1000 X 4 + ^1001 X 4 X 1 + ^1010 X 4 X 2 + /3 1011 X 4 X 2 X 1 
+ %00 X 4 X 3 + ^1101 X 4 X 3 X 1 + %110 X 4 X 3 X 2 + %11 X 4 X 3 X 2 X 1 


In general, a /3 whose subscript in binary notation has ones in the ij, ig, • • • , i^ 
locations from the right is the coefficient of the X. X: • • -X. interaction. 

H ^2 \ 

The set of all 2 n coefficients or parameters form a group B under the appropriate 
operation. In the alphabetic notation this operation (x) is simply commutative multi- 
plication of the letter subscripts with the exponents reduced modulo 2. In the binary 
notation the operation may also be denoted (x) and defined as 




HI HI 1 ’ 

n n- 1 




■ m. ^ Tc k 
1 n n- 


1* 


•k, “ ^d d 
1 n n- 


r 


• d. 


( 5 ) 


where dj = (kj + m.)(mod 2). Thus £ CBA ® /3 DCB = P^^ = ^DA’ and %11 ® 

^lllO = ^loor The definin S Parameter groups that define the fractional replicates at 
the various stopping points are subgroups of B. The aliased sets of parameters at each 
stopping point are the cosets of B(h), which will be denoted Pi ® B(h). 

The principal reason for introducing these notations is that one major problem of 
finding an optimal design is one of finding an optimal matching of design variables to the 
physical variables of the particular experiment. The physical variables in an experi- 
ment will be denoted as Xj, Xg, • • • , X n . It is assumed that the experimenter de- 
cides that these are the only independent variables to be investigated. Each Xj repre- 
sents one of the physical variables and is fixed for the remainder of the experiment. 

For example, 
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Xj = Temperature 
Xg = Time 

X fi = Velocity 

The design variables will be denoted as X^, Xg, X^, . . . and so forth. These 
variables represent abstractions, and tables exist which tabulate experimental designs 
in terms of these design variables. When an experimenter consults one of these tables 
and chooses a design, he must then determine a matching of the design variables and 
the physical variables. Ordinarily the choice is arbitrary because the experimenter 
usually does not have prior information available which would indicate that one matching 
might be preferred to another. A combination of choices of d.p.g. *s, physical- design 
variable matching, and parameter- estimator matching completely specifies for the ex- 
perimenter how to proceed with his experiment and estimation of parameters. Hence, 
such a combination of choices will be called a DESIGN. 

Sidik and Holms (ref. 1) present an analysis of choosing a best DESIGN under the 
following conditions: 

(1) For each 0. of equation (3), the experimenter can specify the probability that 
/3. is not equal to zero, Pj = Pr^ * 0). 

(2) For each h denoting a possible stopping point of the experiment, the experi- 
menter can specify the probability of stopping exactly at the h stopping point, p^. 

(3) For each p^ of equation (3) and each h, the experimenter can specify the value 
to him of obtaining an unbiased estimate of /3^. This is denoted by u^h). 

None of the p's may be separately estimated from a fractional factorial experiment 
unless some assumptions about certain of the P's are introduced. Conditions (1) and 
(3) provide assumptions that will enable the experimenter to assign the estimator for an 
alias set to a single parameter from the alias set and evaluate the consequences of this. 

Changing the matching of physical and design variables will usually change the alias 
sets. For example, if the matching for n = 4 is 
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I 


then the alias set (/3 A , /3g, ^DBA’ * s ma PP e( ^ bito (^0001’ ^0010’ ^1011’ ^1000^ _ 
(/3j, /3 2 , /Sjj, j3g). But the matching 


x i=x B 


X 2= X A 


X 3 =X D 


X 4 = X C 


maps (/3 a , P B , P dba , P d ) into (/3 q 010 , /3qoo 1’ %L11> ^OlOtP “ ^2’ ^ 

Before considering how best to match physical and design variables, let us assume 
that some such matching has been made. Then consider the problem of matching esti- 

j-l 

mators and parameters at the h tn stopping point. The d.p.g. is B(h) and the alias sets 


are 


all those distinct cosets of the form P i (x) B(h) = j/3j , jSj , • • • , /3j 1. 

If the parameter B(h) and the estimator for that alias set is assigned to /3^, 

then, assuming independence, the prior probability that the estimator will be unbiased is 


jei®B(h) 

j*k 


(i - pj) 


th 


Since u^h) is the utility of an unbiased estimate of /J. at the h stopping point 


U(i,k) =u k (h) 


jei®B(h) 
j*k 


o - Pj ) 


(6) 


is the expected utility of the decision to assign the estimator for the alias set p^ ® B(h) 
to the parameter p^. Thus the Bayes strategy is to assign the estimator to the param- 
eter of the alias set which maximizes this expected utility. One case deserves special 
mention. 

Suppose an estimator is confounded with a block effect. It may be safely assumed 
that an unbiased estimate of a block effect has no utility to the experimenter. Let the 
prior probability of the block effect being nonzero be denoted by p^. Then this informa- 
tion can be incorporated into the decision procedure by computing the expected utility as 
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U(i,k) = u k (h) 


jei0B(h) 

j*k 


CL - Pj ) 


(i - P b ) 


(7) 


where = Pr (the block effect does not equal 0). 

With respect to block effects, it is important to note that, depending upon how the 
block parameters are defined, the estimator for the d. p. g. may also be confounded with 
blocks. Let U(i,k ) = max[U(i,k):kei ® B(h)], where the U(i,k) are computed as 

IflaX 

in equation (6) or equation (7), as appropriate. Then, for the assumed physical- design 
variable matching and the given d. p. g. , the maximized expected utility at the n n stop- 
ping point may be denoted by 


u <« =E <«> 

where the summation is over all the distinct cosets at the h™ stopping point. By condi- 
tion (2) (p. 7) it is also assumed that the experimenter can specify the probabilities of 
stopping exactly at each of the stopping points. Thus, 


u = E Psh u < h > <« 

h 

represents the maximized expected utility of the resulting DESIGN. 

The Bayes procedure for choosing an optimal DESIGN is to compute the expected 
utility for each choice of DESIGN and then use any one which yields a maximum expected 
utility. This can be done by computing U as defined in equation (9) for each choice of 
physical- design variable matching and all possible distinct (that is, not equivalent under 
a permutation of letters) choices of d.p. g. 's. NAMER computes U for all possible 
physical- design variable matchings for a specified set of d.p. g. 's. Repeated application 
of NAMER to different choices of d. p. g. 's would then allow the experimenter to carry 
out the full Bayes procedure if he wished. Thus the computing problem is that of mech- 
anizing the evaluation of U for all the matchings of design variables to physical vari- 
ables. 
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ALGORITHM 


The generation and evaluation of all the matchings of the physical and design vari- 
ables present two computing problems. The first problem is the generation of all the 
matchings. This really amounts to computing all permutations of the design variables. 
The second problem is that of evaluating any given permutation. 

The information necessary to evaluate a particular matching is (1) what parameters 
are in the alias sets, (2) the prior probabilities of the parameters, (3) the utility of un- 
biased estimates of each parameter, (4) the alias sets confounded with blocks, and 
(5) the prior probabilities of each block parameter not being zero. The computing pro- 
cedure used by NAMER uses the group properties of the parameters and binary notation 
for the subscripts of the parameters. The parameters arranged in the standard order 
are uniquely identified by the standard- order subscript. Thus, arrays called PROB and 
UTIL may be set up such that the J th entries are Pj_ ^ and Uj_ respectively. Also, 
arrays called BLOCK and PBLOCK may be set up which indicate the alias sets con- 
founded with blocks and the probabilities associated with them. Then when information 
about /3j is needed to compute the expected utilities of equation (6) or (7), it can be 
immediately retrieved. 

For the remainder of this discussion consider only a single stopping point since the 
following procedure will simply be repeated for each stopping point: To determine the 
alias sets, the d.p.g. must be known. Suppose the d.p. g. is stored in an array called 
DPG. The numbers in the array DPG are the standard- order subscripts of the param- 
eters in the d. p. g. when the standard order is computed with respect to the design vari- 
ables. For example, suppose the d.p.g. for the stopping point under consideration is 
j/3j, 0 CBA , /3 dcb , /3 da |, and the matching to be evaluated is 

X 1 = X A 
X 2= X C 
X 3 “ X B 

X 4=X D 

Then the set of standard- order subscripts would be { 0000, 0111, 1110, 1001 } or {0, 

7, 14, 9 }. Since the d.p.g. must always contain the identity or /3j, this is redundant 
information to store. Hence, the numbers which should be stored in DPG are 7, 14, 
and 9. 
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The operation (x) defined by equation (5) can be defined by the Exclusive Or (IEXOR) 
function (defined on p. 24), which is available in almost all computing languages. Thus 
to identify the alias set corresponding to any specified parameter, say /3j, all that is 
needed is to compute the Exclusive Or between J and each number in DPG. The result 
will be the standard- order subscripts of the parameters aliased with /3j. 

If we specify that the numbers in BLOCK are the standard- order subscripts (with 
respect to the design variables) of one parameter from each of the alias sets confounded 
with blocks and that the respective elements of PBLOCK are the prior probabilities of 
the block parameters, then the same use of Exclusive Or can be applied. For example, 
suppose the d.p.g. under consideration is that given previously, j/3j, /^dcB’ 

Pda}- Also suppose it is known that {/3 BA , 0 C , /3 DCA , /3 DB ] and |j 3 dc , 0 DBA , j 3 b , 
(SqA are each confounded with block parameters with prior probabilities of 0. 50. Then 
one element from each of the two alias sets may be chosen to represent it. Suppose they 
are /3 BA and /S B £. Then BLOCK(l) should be set to 11) 2 = 3, and BLOCK(2) should be 
set to 1100) 2 = 12; and PBLOCK(l) = PBLOCK(2) = 0. 50. 

It is now an easy task to compute the expected utilities of equations (6) and (7). It 
remains to find all the distinct alias sets in some economical manner. To do so, set up 
an array denoted T1 which is 2 n words long. This will be used as an indicator array to 
indicate if a parameter has been found in an alias set so far. Begin the computation for 
the stopping point by setting U(h) = 0. 0 and initializing the T1 array to some value, say 
zero. 

For each block effect set the element given by BLOCK in T1 to some indicator value 
not equal to the initialization value, say IRUN; and compute the Exclusive Or of that 
element and every value in DPG. This will yield the standard- order subscripts of all 
the parameters in the alias set. To indicate that these parameters have been identified 
as members of an alias set, set the locations of T1 corresponding to these parameters 
equal to IRUN. These standard- order subscripts and the value in PBLOCK indicate 
where to find the probabilities and utilities necessary for making the optimal estimator- 
parameter matching according to equation (7). Compute the expected utilities, identify 
the maximum, and add this value to U(h). 

Now begin searching T1 until a value not equal to IRUN is found. Suppose the sub- 
script of this value is K. Then compute the EXOR of K - 1 with each number in DPG. 
Along with K - 1 itself, this will yield the standard- order subscripts of all parameters 
in the alias set containing To indicate that these parameters have been identified 

as members of an alias set, set the locations of T1 corresponding to these parameters 
equal to IRUN. 

These standard- order subscripts provide the information needed to find the proba- 
bilities and utilities necessary for making the optimal estimator-parameter matching. 
Compute the expected utilities, identify the maximum, and add this value to U(h). 
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Now continue searching T1 from the location K + 1 for another value not equal to 
IRUN. This will find the next parameter in the standard order which has not yet ap- 
peared in an alias set. Thus, the preceding evaluation procedure should be repeated 
until the end of the T1 array is reached. At that point, all the distinct alias sets will 
have been identified and evaluated once and only once. The value of U(h) will then be the 
total maximized expected utility for the h th stopping point corresponding to the optimal 
estimator- parameter matchings for the current physical- design variable matching. 


This same procedure is simply repeated for each stopping point and then 


u =Z e S h u < h > 


may be calculated. The matchings of physical- design variables which provide the lar- 
gest values of U(h) and U may be kept updated in several arrays. Then when all the 
permutations are completed, the optimal matchings will be available. 

What must now be developed is a procedure for generating all the matchings of phys- 
ical to design variables in some economical manner. Since all the permutations are to 
be evaluated, the result does not depend upon which matching is done first or in what 
order they are generated. Thus the starting permutation and the order of generation 
may be whatever is most convenient computationally. A simple convention used in 
NAMER is to begin with the matching 


X 


1 


= x, 


X 


2 



The distinction has been made previously that the alphabetic subscripts are for the de- 
sign variables and the numeric subscripts for the physical variables. Thus the preced- 
ing starting convention has both sets of variables in the standard ordering. Suppose the 
d. p.g. at a given stopping point is |/3j, Pq-qj^, ^DCB’ ^Da}' Then for the matching 


X 1= X A 


X 2 =X B 


X 3 =X C 


X 4= X D 


( 10 ) 


the DPG array contains 0111 for (Pq#a)’ 1110 for (j3 D ^ B ), and 1001 for (/3 DA ). To 
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evaluate a different matching, say 


X 1= X A 


X 


2 



X 3= X D 


X 4 = X B 



(U) 


the DPG array should contain 1011 for (Pq-qa)> 1110 for (/Sg^-g), and 0101 for (/3g A )- 
The latter DPG can be derived from the former by permuting the binary digits according 
to the same permutation that gives the ordering X B , Xg, Xq, X^ starting with Xg, X^, 
Xg, X A . Recall that the binary digits are numbered from right to left. Thus the dif- 
ferent matchings of variables can be achieved by constructing all the n! permutations 
of the rightmost n binary bits in the numbers in DPG. The same procedure applies to 
the BLOCK array for the same reasons. 

Or d- Smith (ref. 4) has presented a survey of a number of possible permutation 
algorithms. Of these, the best for the purposes of NAMER is the one by Trotter 
(ref. 5). Trotters' algorithm computes all the permutations as a sequence of adjacent 
transpositions. To see why this is best, consider how to achieve the permutation of the 
binary bits by means of arithmetic and logical machine operations. Let M be the num- 
ber to be changed and express it in binary as M = m n m n _i- . . m^m^. . .mj and 
suppose the digits mj and nijj are to be transposed. Compute 


n 



j 

j-1 

1 


I = AND (0 

0 

0 . 

. . 1 

0 . . 

. 0, 

M) 

n 



j 

j-1 

1 


J = AND (0 

0 

0 . 

. . 0 

1 . . 

. 0, 

M) 

n 



j 

j-1 



K = AND(1 

1 

1 . 

. . 0 

0 . . 

. 1, 

M) 


I = 1/2 


J — 2 

M=I + J + K = in in n . . . m- ,m. . . . m, 
n n-1 ]-l ] 1 
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Notice that the shifting of the digits nij and mj_ ^ is accomplished by the multipli- 
cation and division by 2. If the permutations were not the result of transpositions of 
adjacent digits, a more general shift function would be needed or the use of powers of 2 
would be needed. These would take more time to execute and/or more logic to control 
than the current method. This is an important consideration since the computing of the 
permutations accounts for a substantial portion of the computing job. 

A third major problem involved in the program is that of providing the necessary 
output from the calculations in an economical and useful manner. To explain what 
NAMER does it will be convenient to introduce some notations for, and properties of, 
permutations. Let A denote the set of the first n letters of the alphabet arranged in 
order; that is, A = { A, B, C, D, . . . }. Let an ordering of A be the set of the first 
n letters of the alphabet arranged in some arbitrary order. Let a permutation on A or 
any particular ordering of A be a function denoted as 

/I 2 3 . . . n \ 

P = (12) 

V 1*2*3 • * • V 

Xl_ 11 

which means to take the j element of the ordering and make it the ij element for 
j = 1, 2, . . . n. Thus a permutation is a function which maps the set of all possible 
orderings of A one-to-one onto itself. Since the upper line of equation (12) is redun- 
dant, the notation for P is often reduced to P = (ij, ig, . . ., i R ). 

A transposition is a permutation which interchanges exactly two elements of A. 

Any permutation can be expressed as a product of transpositions of the form (1, ij) 

(2, ig) . . . (n, i n ). Here (j, i- ) is an abbreviated notation for 

1 2 . . . j . . . ij . . . n 

1 2 . . . Ij . . . j . . . n 

The product of transpositions may be expressed as a transposition vector (ij, . . . , i n ) . 
As an illustration note 
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= (1, 3)(2, 3)(3, 5)(4, 5)(5, 5) =<3, 3, 5, 5, 5> 

Then P[ABCDE] = [BDAEC] directly and 

(1, 3)(2, 3)(3 , 5)(4, 5)(5, 5)[ABCDE] =(1, 3)(2, 3)(3, 5)(4, 5)[ABCDEl 

= (1, 3)(2, 3)(3, 5)[ABCED] 

= (1, 3)(2, 3)[ABDEC] 

= (1, 3)[ADBEC] 

= [BDAEC] 

as a sequence of transpositions. This illustrates the equivalence of the two ways of ex- 
pressing P. It is obvious that the permutation expressed as a product of transpositions 
is most convenient for this computer application. Let NDPG be the number of defining 
parameter groups. Then the orderings O^, . . . , OnDPG +1 which yield the largest 
overall expected utility (Oj) and the largest utilities at each stopping point (Og, . . . , 
C>ndpg+i) are the information the statistician seeks. 

The program included in this report does not print all the information for every per- 
mutation since this would be much too large a volume of output. All that is saved are 
the orderings O. (i = 1, 2, . . . , NDPG + 1). The initial letter- factor matching is the 
assignment of the i in letter of the alphabet to i. After all the permutations are per- 
formed and evaluated, the program returns the order of the letters and the d. p. g. *s to 
their original state. 
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Then permutations P^^ must be found to effect reorderings upon A = {A, B, C, D, 
E, ... } to achieve Oj, 0 2 , . . . °ndCG+ 1 in some efficient manner. Let 

P t [A] = Oj 

p 2 I a 1 = = °2 

P i[A] = Pi^rAtFi-ltAD] = °i 

and define O q = A, P Q = (1, 2, . . M). Then the sequence of permutations 

PjJpT^jJO.-jJ should be determined. If P = (p-p . . p m ), then P" 1 = (r^, . . 

r ), where r = i. That is, r, is the subscript of the location in P which contains j. 
m Pi 1 J , 

If Q = (q 1 , . . q m )and P _i = (rp . . r m ), then QP = (Sp . . s m ), where 

s i = q r . That is, to get Sj, find the character in the j th position of P” and then go 
to that location in Q to find Sj . Putting the two operations together - if 

P = (Pp • • • > Pjjj) 

Q = (qp • • • , q m ) 

then 

QP' 1 = (s x , . • • , s m ) 

where Sj = q r . That is, find the subscript of the location in P which contains i and 
go to that location of Q to find s^. 

PROGRAM DESCRIPTION 

The NAMER program is composed of the main program NAMER and five subrou- 
tines: BLOK, LINER, PERMUT, REMACH, and RECT. They will be discussed in 
some detail after the following brief descriptions: 
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Program name Calling name 


Purpose 


NAMER - 

Main program 

Input; evaluation of each permutation; identifies and 
and saves optimal matchings; overall program 
control. 

BLOK 


Block data subprogram. 

PERM 

PERMUT 

Determines permutations and permutes DPG and 
BLOCK arrays. 

REMAX 

RE MACH 

Achieves rematchings of physical- design variables; 


outputs detailed description of matchings of 
physical- design variables and estimators- 
parameters by appropriate calls to LINE and 
ALINE. 


LINER LINE 


ALINE 


ERECT RECT 


Prints one line of output identifying an estimable 
parameter and the utility of the assignment of the 
estimator to the parameter. 

If two or more members of an alias set are tied for 
maximum utility, ALINE identifies the remainder 
not printed by LINE. 

Prints summary output table of the Bayes matching 
of physical to design variables, the optimal match- 
ings for each stopping point, and the utilities. 


Main Program NAMER 

NAMER is the main program and is divided into 16 major sections, as indicated by 
the comments cards in the listing. 

Sections 1 to 10 . - Read and write input information. 

Section 11 . - This section is the heart of the program, where each d.p. g. is eval- 
uated and its contribution to the total expected utility computed for a given permutation 
of the letters. The section is divided into two subsections, 11A and 11B. Section 11A 
chooses the parameter-estimator matching for the alias sets which are confounded with 
blocks. Section 11B does the same for the remainder of the alias sets. The computa- 
tions for IUTILF = 1, 2 are less complicated than those for IUTILF = 3, 4, 5. Thus 
these cases are separated in the program. See section 7 of the input description for 
further explanation of IUTILF. 


17 


Section 12 . - The expected utility of this ordering (PBAYEX) and the expected utili- 
ties at the stopping points (PSUMX(- )) are compared to the best utilities to this point 
(PBAYES AND PSUM(-)). If any one or more is larger than the best so far, the current 
appropriate utility is placed in PBAYES or PSUM(-) and the ordering (indicated by the 
contents of IALPHA(-)) is saved in ISAVEP(-, I). The convention is that ISAVEP(-, 1) 
saves the ordering which gave the best weighted utility and ISAVEP(-, I + 1) saves the 
ordering which gave the best utility for the I m d. p. g. 

Section 13 . - This section permutes the current letter- variable relation for the 
appropriate class or classes. If all possible distinct permutations have been realized, 
control passes to section 16. If not all permutations have been realized, control passes 
to section 14. 

Section 14 . - If current execution time exceeds that allowed, go to section 15. 
Otherwise go to section 11 to evaluate the current ordering. 

Section 15 . - All essential information is punched on cards to permit a restart of 
this case on another computer run beginning with the current ordering. 

Section 16. - Print current clock time and call REMACH. 


Subroutine PERM 

This subroutine uses a FORTRAN translation of Trotter’s routine (ref. 5) to per- 
mute the elements in an array as a sequence of transpositions of adjacent elements. The 
first two blocks are the logic which determines which two adjacent elements are to be 
transposed. The third block (beginning with line 44) is where the arrays IALPHA, DPG, 
and BLOCK are actually permuted. 


Subroutine REMAX 

This subroutine uses the information saved in section 12 of NAMER to recompute 
the utilities of the optimal matchings. 

Section 1 . - At this point, the DPG and BLOCK arrays contain the same values they 
had as initial input data. The permutation required to achieve the first optimal ordering 
is computed and stored in KPERM(-). 

Section 2 . - Write the letter- variable matching. 

Sections 3 and 4 . - Translate permutation vector into transposition vector, and per- 
mute BLOCK and DPG arrays. 

Section 5 . - Performs same function as section 11 of NAMER except that calls to 
LINE and ALINE are made as appropriate. 
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Section 6 . - Output overall expected utility and expected utilities for each d. p. g. 
Section 7 . - Compute next permutation and shift to section 2. 

To illustrate sections 1, 3, and 7 of REMAX consider the following example: 

Oj - [4 1 3 2 5] = Pj[A] 

0 2 = [3 4 1 5 2] = P 2 [A] 

Og = [4 1 3 2 5] = P 3 [A] 

0 4 = [4 3 2 1 5] = P 4 [A] 

O s = [1 4 5 2 3] = P 5 [A] 

The sequence of values is as follows: 


Sequence 

KPERM 

KKSAVE 

K2CYCL 

Ordering 

1 

P x = (24315) 

P x = (24315) 

{ 24345 } 

[41325] 

2 

P 2 = (35124) 


(23154) = PjjPj 1 


3 

PgPj 1 = (23154) 

P 2 = (35124) 

{ 23355 ) 

[34152] 

4 

P 3 = (24315) 


(31254) = PgP^ 1 


5 

PgPj 1 = (31254) 

P 3 = (24315) 

{ 33355 ) 

[41325] 

6 

P 4 = (43215) 


(14235) = P 4 P 3 X 


7 

P^’ 1 -- (14235) 

P 4 = (43215) 

{ 14445 ) 

|43215] 

8 

P 5 = (14523) 


(25413) = PgP^ 1 


9 

PgP^ 1 = (25413) 

P 5 = (14523) 

{ 25455 ) 

(14523| 


Subroutine LINER 

This is a double- entry subroutine with entry points LINE and ALINE. LINE is 
called by output once for each alias set. The entries of the calling vector are the 
standard- order subscript number of the parameter in the alias set to which the esti- 
mator should be assigned and the expected utility of that assignment. The standard- 
order subscript is then used to identify the parameter in terms of the interaction of the 
independent variables it measures. This identification is then printed in numerical form 
and in Hollerith form using the first six characters of each factor identification card. 

The utility is also printed. 
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ALINE is called whenever there are two or more parameters aliased which give the 
same maximized expected utility. The call to LINE causes one of the aliased param- 
eters to be identified and the expected utility of the estimator to be printed. The call to 
ALINE causes the remaining parameters to be identified. They are, however, identi- 
fied only by the numerical form of their interaction. 


INPUT DESCRIPTION 

The following is a detailed description of the input data necessary to run a problem. 
There are nine basic sets of input data. Each is described in detail here. An example 
of the type of problem to which this program may be applied is given in reference 1 and 
a similar problem is discussed in appendix A. A sample set of data for this problem is 
given in table I. A pictorial illustration of the input deck setup is given in figure 1. 
Multiple cases may be run back-to-back. The last card of the last case should have 
ENDALL punched in the first six columns. 

The nine basic sets of input data are as follows: 

(1) IDENTIFICATION (13A6, A2) (IDENT). This is one card; all 80 columns are 
used for Hollerith identification of problem. 

(2) MAXIMUM TIME (F6. 0) (TMAXX). This is the maximum machine time in min- 
utes permitted for this case. If this time is exceeded and the case is not fully evaluated, 
all pertinent information is punched on cards to permit a restart of the program. 

(3) TYPE OF RUN (16) (ITYPRN). 

(3A) SPECIFIED MATCHING (9A1) (XT3). A "1” for ITYPRN indicates this is a 
regular first-time rim and data sets 4 to 9 will be read. A "2" indicates this is a re- 
started case and only those cards punched by the previous run need be read. A "3 M 
indicates that only one matching will be evaluated. This matching is specified on the 3A 
card with the first n letters of the alphabet (excluding I) in the first n columns of the 
card in the order appropriate to the matching desired. 

For example, if to indicate an evaluation of the following matching is desired 

X 1= X B 
X 2= X A 
X 3= X D 
X 4 = X C 

one card would be supplied with BADC in the first four columns. 
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(4) NUMBER OF FACTORS (16) (NFAC). Up to nine factors can be considered. 

(5) FACTOR IDENTIFICATIONS (13A6, A2) (FAC). One card for each factor. The 
first six characters of each card are used as output identification, so they should serve 
as useful abbreviations. 

(6) NUMBER OF CLASSES (16) (NCLASS). 

(6A) NUMBER IN EACH CLASS (916) (NSUBI). NCLASS is the number of classes of 
factors. If this is 1, input set 6A is not read. If the number of classes is more than 
one, card 6A specifies the number of factors in each class. The factors within a class 
will be permuted among each other, but permutations between classes will not be per- 
mitted. The first NSUBI(l) factors will be assumed to belong to the first class, the next 
NSUBI(2) to the second class, and so forth. Holms and Sidik (ref. 6) present an experi- 
ment in which there are two classes of variables which could not be mixed. Most ex- 
periments have only one class. 

(7) NUMBER OF NONZERO PROBABILITIES, UTILITY FUNCTION, AND CON- 
STANT (NPIN, IUTILF, UCOEF) (216, F10. 9). The number of parameters with nonzero 
prior probabilities is specified in the first six columns. The choice of utility function is 
indicated in the second six columns. UCOEF is used in defining utility function 5 and is 
given in the next 10 columns. Each parameter with a nonzero prior probability or utility 
is identified in terms of the integer subscripts of the independent variables in the inter- 
action with which it is associated. See the input set (8) description for further informa- 
tion. The possible choices of utility function here are 

(a) IUTILF = 1 

{ 1. 0 for unbiased estimators 
0. 0 for biased estimators 

(b) IUTILF = 2 


u i = 


'p^ for unbiased estimators 
0. 0 for biased estimators 


(c) IUTILF = 3 


u. 

l 


{ Xj for unbiased estimators 
0. 0 for biased estimators 


where Xj is given with the p^ in input set (8). 
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II III 


III II II 


(d) IUTILF = 4 

f p^x^ for unbiased estimators 
1 [ 0. 0 for biased estimators 

where Xj and Pj are given in input set (8). 

(e) IUTILF - 5 

{ UCOEF • Xj + (1 - UCOEF)p i for unbiased estimators 
0. 0 for biased estimators 

where and p. are given in input set (8) and it is assumed 0. 0 ^ UCOEF ==1.0. 
It should be noticed that these utility functions do not depend upon the stopping point 
as implied by condition (3) on page 7. To provide a capability of making the utility func- 
tion depend upon the stopping point, the user can weight the stopping points by use of the 
weighting values WT(I) read in input set (9B). Thus the m(h) used in the program are 
computed as u^(h) = u^ * WT(h). 

(8) PRIOR PROBABILITIES AND UTILITIES (911, 2F10. 0) (IT1, P, UT). Each 
parameter with nonzero prior probability is identified in terms of the integer subscripts 
of the independent variables in the interaction with which it is associated. These sub- 
scripts may be supplied in any order anywhere in the first nine columns of the card. 

The prior probability and the utility follow with 10 columns each, in F10. 0 format. The 
utility need not be specified if IUTILF = 1 or 2 as previously described, for the program 
then supplies the utility. If IUTILF = 3, 4, or 5, the utility must be specified explicitly. 
For example, suppose the XgXgX^ interaction parameter = /3 22 is assumed to 

satisfy P(0 22 ^ 0) = 0. 850 with utility of 0. 95. Then the card input could be 
bbb5b2b3b. 850bbbbbb. 95. If the prior probability of a parameter being nonzero is zero, 
no data need be supplied for that parameter. 

(9) NUMBER OF DEFINING GROUPS (16) (NDPG). For each d.p. g. (as many as 32 
permitted) there must be one set of inputs 9A to 9E. 

(9A) IDENTIFICATION OF STOPPING POINT (4A6) (IDDCG). 

(9B) NUMBER OF GENERATORS, PRIOR PROBABILITY OF STOPPING, WEIGHT- 
ING VALUE (16, F6.6, F6. 0) (NGEN, PSTOP, WT). If the d. p. g. corresponds to a 
(l/2) r fractional replicate, r independent generators must be supplied. Program limi- 
tations restrict NGEN to values less than or equal to seven. 

(9C) THE GENERATORS OF THE d.p.g. (9A1). The generators are supplied in 
terms of the first NFAC letters of the alphabet on the first nine columns of the card. 
There is one card per generator. For example, if the d.p.g. at a particular stopping 
point is {/3 p /3 edcba , /3 cba , /3 ed , 0 da , 0 ECB , 0 DCB , 0 EA }, three generators are 
sufficient; and one such choice might be /3^ bA , %£>> and /3 DA . Three cards which 
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will define the above d. p. g. might be 


AbBbC 

ED 

bbAbbD 

The order and position of the letters is unimportant as long as they are on the first nine 
columns of the card. If the number of generators is zero, no type-9C cards are read. 

(9D) NUMBER OF BLOCK PARAMETERS (16) (NBLOCK). 

(9E) IDENTIFICATION OF ALIAS SETS CONFOUNDED WITH BLOCK EFFECTS 
AND THE PRIOR PROBABILITY ASSOCIATED WITH THE BLOCK EFFECTS (9A1, 

F5. 5) (XT3, PBLOCK). Any single parameter from an alias set which is confounded 
with a block effect may be input in terms of the first NFAC letters of the alphabet in the 
first nine columns of one card. This is followed by the prior probability of the block 
effect on the next five columns. There is one card for each block effect that has a non- 
zero prior probability. If the alias set is the d. p. g. , the first nine columns may be left 
blank. 


OUTPUT DESCRIPTION 

The first part of NAMER output is the printout of the input data. This is followed by 
NDPG + 1 printouts. The first set is for the Bayes DESIGN which optimizes the overall 
expected utility. The subsequent sets are for the DESIGNS that optimize the expected 
utilities for the individual stopping points. Each of these sets of output consists of the 
following: 

(1) The optimal matching 

(2) Tables of the parameters chosen to be estimated and the expected utilities of 

these choices 

(3) The overall expected utility and the expected utilities at the stopping points 
For example, the first two pages of the sample output in appendix A indicate the 

input data. The next two pages provide the information about the Bayes DESIGN, as the 
label indicates. The Bayes matching is seen to be 

X^TEMP) = X c 

x 2 (press) = x D 
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X 3 (TIME) = X B 


X 4 (VEL) =X E 
X 5 (ANGLE) - X A 

Then for d. p. g. number one (which is the 1/ 4 replicate of the full factorial) the choices 
of parameters to be estimated are indicated. Each parameter is identified in terms of 
the integer subscripts of the independent variables in the interaction with which it is 
associated. They are also identified by the Hollerith identifications input in section (5) 
of the input, and the utility of the choice is printed at the far right of each line. Thus 
the first line of output for d. p. g. number one indicates that the coefficient of 
X ^ 2 X 3 = X C X D X B has been chosen from its alias set as the parameter to be estimated. 
This interaction is the TEMPxPRESSxTIME interaction, and the expected utility of 
this choice is 0. 20. This utility value does not include the weighting factor at this point. 

Below the detailed output for the three d. p. g. 's are printed the overall expected 
utility and the utilities for each of the d. p. g. 's for this matching. 

Similar output provides the detailed output for the designs which maximize the ex- 
pected utilities for each of the stopping points. The format and arrangement are the 
same as that described for the Bayes matching. This is followed on the last page by a 
summary table providing the various matchings, their expected overall utilities, and 
expected utilities at the stopping points. 


SPECIAL LEWIS RESEARCH CENTER ROUTINES 

Some of the following functions and subroutines available in the FORTRAN IV - 
Version 13 language at the Lewis Research Center may not be available (or not available 
in FORTRAN) at other computer installations. Thus, their usage is explained, and the 
user can write functions or subroutines providing the same capabilities in a language 
compatible with the available computer. 

The functions and subroutines available at Lewis are the following: 

(1) AND(A, B). A real function of the Real or Integer variables A and B. Like 
bit positions of A and B are compared. A 1 is placed in those positions of the result 
where there are l's in both A and B, and a zero is placed in the result otherwise. 

(2) IEXOR(A, B). An integer function of the Real or Integer arguments A and B. 
Like bit positions of A and B are compared. A 1 is placed in those positions of the 
result where exactly one of A or B is a 1, and a zero is placed in the result other- 
wise. 
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(3) IALS(N, X). An integer function of Integer N and Real or Integer X. The con- 
tents of X are shifted to the left N places and zeros put into the vacated rightmost 
positions. 

(4) IARS(N, X). An integer function of Integer N and Real or Integer X. The con- 
tents of X are shifted to the right N places and zeros put into the vacated leftmost 
positions . 

(5) BCREAD(X1, X2) and BCDUMP(X1, X2, K). These subprograms provide for 
input and output in absolute binary. A call to BCREAD(X1, X2) causes cards to be read 
in binary format at the rate of 22 words per card. The data are stored sequentially in 
the core, beginning with the address of the variable XI and ending with the address of 
the variable X2. A call to BCDUMP(X1, X2, K) causes cards to be punched in binary 
format at the rate of 22 words per card. The data are taken sequentially from the core, 
beginning with the address of variable XI and ending with the address of variable X2. 

K provides card numbering control and is always set to zero by NAMER. 

As an example of the usage of these routines, consider the first call to BCREAD in 
section 10 of NAMER. DUMPl(l) is equivalenced to NFAC. NFAC is the first variable 
of nine variables in the labeled common block Bl. LD1 is set to 9 at the start of 
NAMER. Thus, the call BCREAD (DUMPl(l), DUMPl(LDl)) causes the variables 
NFAC, NCLASS, NN, NDPG, PBAYES, and so forth, to be read from unit 5 in binary 
format. 

(6) TIMEl(X). This subroutine enables the programmer to read the storage cell 
clock. The following illustrates the procedure for using TIME1 to calculate elapsed 
time. 


CALL TIME 1 (XI) 


CALL TIME1(X2) 


Then 


X2 - XI = Clock pulses 


X2 - XI 

= Elapsed time in seconds 

60 


X2 - XI 

= Elapsed time in minutes 

3600 
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(7) OR (A, B). A real function of the real or integer arguments A and B. Like 
bit positions of A and B are compared. A 1 is placed in those positions of the result 
where either or both of A and B are a 1. A zero is placed in those positions of the 
result wherever both A and B are zero. 

TIMING INFORMATION 

Several sample problems using single telescoping were run on NAMER to estimate 
the amount of time required by the program. For these problems, the first stage was 
assumed to be the smallest experiment large enough to estimate all main effects, and 
the last stage was the full factorial. Each problem was run once using utility function 2 
and once using utility function 3. The results are summarized in table II and figure 2. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 8, 1971, 

132-80. 
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APPENDIX A 

SAMPLE PROBLEM AND PROGRAM OUTPUT 

Consider a five- factor experiment involving 

= Temperature 
Xg = Pressure 
Xg = Time 
X 4 = Velocity 
X & = Angle 

Suppose that the experimenter's facilities are such that he can only perform four treat- 
ment combinations at one time and be reasonably sure that experimental conditions are 
homogeneous. Thus his experiment should be designed as a blocked factorial design 
with blocks of size four. Assume also that he has enough materials at one time to per- 
form eight treatment combinations, but no more, and that batches of uniform material 
are not available in quantities that will supply more than eight treatment combinations. 
Then the blocks of the experiment might be as shown in the illustration, where the two 
columns represent two different test facilities and the four rows represent four different 
batches of raw material. 


Column blocks 


Row blocks 
(batches) 


(test facilities) 


^ 1 - 

2 ' 

1,1 


I 1 ’ 2 ! 



2,1 


Rl 

CSl 



3,1 


| 3,2 | 



I 4 ’ 1 


4,2 j 
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The difference between the first block and the second block in a row is due to perform- 
ing the experiment in two different test facilities. The differences between rows are 
due to possible effects of new batches of materials. Suppose the experimenter feels that 
there is a probability of 0. 50 of there actually being a test facility block effect. Let the 
probability of there being an effect due to differing batches of raw materials be 1.0. 
Assume further that the probability of an interaction between these block effects is spe- 
cified as zero. The stopping points of the experiment at each stage are 

(1) Stage one, after completion of blocks (1, 1), (1, 2) 

(2) Stage two, after completion of blocks (1, 1), (1,2), (2,1), (2,2) 

(3) Stage three, after completion of the full factorial 

Based upon his available resources and upon past histories of some similar projects he 
has worked on, the experimenter feels probabilities in the following table are appro- 
priate: 


Coefficient of- 

Standard- order 

Prior probability of 


subscript 

being nonzero 

x o 

0 

1 . 

00 

X 1 

1 


80 

X 2 

2 



X 2 X 1 

3 



X 3 

4 



X 3 X 1 

5 



X 3 X 2 

6 



X 3 X 2 X 1 

7 



X 4 

8 

1 . 

0 

X 4 X 1 

9 

.50 

X 4 X 3 

12 

.50 

X 4 X 3 X 1 

13 

.40 

X 5 

16 

1 . 

0 

X 5 X 1 

17 

.40 

X 5 X 3 

20 

.30 

Stopping probabilities: p^ g = 0. 30, p 2s = 0. 

P 3s = °' 3° 

40, j 

[ 


All the other coefficients have zero prior probability. 
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Assume the purpose of the experiment is to maximize the response. Also assume 
that the cost of the experiment is about proportional to the number of treatment combin- 
ations rim. Then a reasonable choice for utility function might be 

{ p./n. if unbiased 
1 h 

0 if biased 

where n^ is the number of treatment combinations. To achieve this using NAMER, 
use utility function 2 and for the weighting values at the stages use 1/n^. 

Rather than investigating all the possible nonequivalent d.p. g. 's and their telescop- 
ing options, the best matchings of physical- design variables and parameters to estima- 
tors will be determined for the following choices of d.p.g. 's: 

B(l) (Pp ^CBA’ % DC’ %DBa} 

B(2) ={/ 3 p %DBa} 

B(3) =|/3j| 

Using the d.p.g. |/3 p P CBA , /3 D CB , ^DA’ ^EDC’ ^EDBA’ %B’ ^ECa} for block ^ 
and the rules presented in reference 3, it may be shown that the following assignment of 
treatment combinations will lead to the block confounding presented in the table: 

Block Treatment 

(1, 1) (1), dca, ecb, edba 

(1,2) ba, deb, eca, ed 

(2, 1) db, eba, edc, ea 

(2, 2) da, c, edeba, eb 

(3, 1) a, dc, eeba, edb 
(3,2) b, deba, ec, eda 

(4, 1) dba, cb, edca, e 

(4, 2) d, ca, edeb, eba 
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Stage 

Alias sets confounded with - 

Test facility effect 

Raw material effect 

1 

{^DA’ ^DCB’ @ECA’ %b} 

{/3p 0 CBA , 0 EDC ’ @EDb} 

2 

{%A> %b} 

1 

1 

h’ ^edba} 
^CBA’ %Dc} 

3 

{^da} 

: 

jl 

'M , 

^CBAJ 

/edba} 

[^EDc} 


The sample FORTRAN data sheets given in table II supply the data necessary to run this 
problem as described. The sample output for this problem follows. 


NAMER OUTPUT NAMER SAMPLE PROBLEM 

PROGRAM WILL DUMP FOR RESTART IF NOT FINISHED IN 2. MINUTES. 


CURRENT EXECUTION TIME 0.00 

THERF APF 5 FACTORS. THEY ARE... 


3 

A 

5 


TEMP 

SOURCE 

TEMPERATURE 

PRFSS 

SOURCE 

PRESSURE 

TIME 

TIME DURATION 

VEL 

SOURCE 

VELOCITY 

ANGLE 

ANGLE 

OF INJECTION 


IS PARAMETERS WITH NON-ZERO PRIOR PROBABILITIES AND UTILITIES 


UTILITY FUNCTION 2 


0 

1 . 0 C 0000 

1. 000000 

I 

0.8 00000 

0 . 800000 

2 

0.800000 

0 . 800000 

1 ? 

0.800000 

0 . 800000 

3 

0.800000 

0 . 800000 

13 

0 o 800000 

0.800000 

23 

0. 800000 

0.800000 

i 23 

C . 800000 

0 . 800000 

A 

1 . 0 C 0000 

I . 000000 

H 

0.500000 

0 . 500000 

34 

C . 500000 

0 . 500000 

134 

0.400000 

0.400000 

5 

1.000000 

1 , 000000 

1 S 

0.400000 

0.400000 

35 

0.300000 

0.300000 


? DEFINING PARAMETER GROUPS 


1/4 R FP— — ROW 1 
DP G 1 

? GENERATORS 

PROB OF STOPPING 0.30000 WEIGHT 0.125000 

ABC 
COE 

THERF ARE 2 BLOCK PARAMETERS 

AD 0.50000 

1.00000 
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1/7 REP— ROWS 1,2 


OP G ? 

1 GENERATOR S 

PROB OF STOPPING 0*40000 WEIGHT 0*062500 

ABDE 

THERE ARE 3 BLOCK PARAMETERS 

AO Oo 50000 

To 00000 
ABC I„00000 


FULL— ALL ROWS 
OP G 3 

0 GENERATORS 

PROB OF STOPPING 0.30000 WEIGHT 0*0*312 50 

THERE ARE 5 BLOCK PARAMETERS 

AO C, 50000 

to 00000 
ABC 1.00000 

ABOE 1.00000 

COE l. 00000 

CURRENT EXECUTION TIME 
CURB ENT EXECUTION TIME 


0, 01 
0*03 



OFF I N TNG PARAMETER GROUP NO. 2 
1/2 REP — ROWS l,P 

************************************************************ 


*3 *4 



TIME 

VEL 


C. 250000 

0* 

G KEAN 




c 

*1 *3 *5 

TEMP 


TIME 


ANGLE 

c 

*i 

TEMP 





C.BOOOCO 

*2 


PRESS 




c.eooooo 

*1 *2 

TEMP 

PRESS 




C. 800000 

*3 



TIME 



C. 800000 

*1 *3 

TEMP 


TIME 



C. 800000 

*? *3 


PRESS 

TIME 



C. 800000 

*1 *2 *3 

TEKP 

PRESS 

TIME 



C. 800000 

* 4 




VEL 


1.000000 

*1 *4 

TEMP 



VEL 


C. 500000 

*3 *5 



TIME 


ANGLE 

C. 300000 

*1 *3 *4 

TEMP 


TIME 

VEL 


0.400000 

*5 





ANGLE 

1.000000 

*1 *«j 

TE MP 




ANGLE 

C. 400000 

OEFINING PARAMETER GROUP NO* 3 

FULL— ALL ROWS 

***************** 

* * * * 

* * * 

* * * 

* * * 

»*********< 

K *********** 

*2 *5 


PRESS 



ANGLE 

0 

0* 

G KEAN 




c 

*1 *3 *5 

TE KP 


TIME 


ANGLE 

c 

*2 *3 *4 *5 


PRESS 

TIME 

VEL 

ANGLE 

c 

*1 *? *4 

TE KP 

PRESS 


VEL 


c 

*» 

TEKP 





c.eooooo 

*2 


PRESS 




C. 800000 

*! ♦? 

TEKP 

PRESS 




C. 800000 

*3 



TIME 



C. 800000 

*1 *3 

TEMP 


TIME 



C. 800000 

*2 *3 


PRESS 

TIME 



C. 800000 

*! *2 *3 

TEKP 

PRESS 

TIME 



c.eooooo 
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I 




*4 





VEL 


1.000000 

*1 

*4 


TEMP 



VEL 


C. 500000 

*2 

*4 



PRESS 


VEL 


c 

*3 

*4 




TIME 

VEL 


C. 500000 

*T *3 

*4 


TEMP 


TIME 

VEL 


C.400000 

42 * 3 

♦ 4 



PRESS 

TIME 

VEL 


C 

*1 42 *3 

*4 


TEMP 

PRESS 

TIME 

VEL 


c 


*5 






ANGLE 

1.000000 

*1 



TEMP 




ANGLE 

C.400000 

*1 *2 

*5 


TEMP 

PRESS 



ANGLE 

0 

*3 

*5 




TIME 


ANGLE 

C. 300000 

*2 *3 

*5 



PRESS 

TIME 


ANGLE 

0 

*1 *2 *3 

*5 


TEMP 

PRESS 

TIME 


ANGLE 

c 

• 

*4 *5 





VEL 

ANGLE 

0 

*1 

*4 *5 


TEMP 



VEL 

ANGLE 

c 

*2 

*4 *5 



PRESS 


VEL 

ANGLE 

c 

*1 *2 

*4 *5 


TEMP 

PRESS 


VEL 

ANGLE 

c 

*3 

*4 *5 




TIME 

VEL 

ANGLE 

c 

*1 *3 

*4 *5 


TEMP 


TIMt 

VEL 

ANGLE 

c 

*1 *2 *3 

*4 *5 


TEMP 

PRESS 

TIME 

VEL 

ANGLE 

c 

FOR THE ABOVE PERMUTATION 

THE 

EXPECTED UTILITY 

IS 

0.42169 


THE EXPECTED UTILITIES AT 

THE 

STOPPING POINTS ARE*. 




dffining 

PARAMETER GROUP 

1 

0. 31500 





DEF INING 

PARAMETER GROUP 

2 

0. 59062 





DEFINING 

PARAMETER GROUP 

3 

0.30312 






THIS MATCHING MAXIMIZES THE EXPECTED VALUE AT THE 1 STOPPING POINT 1/4 REP— ROW 1 


VAR! ABLE 

! TFMP 
? PRESS 
3 TIME 
A VFL 
6 ANGLE 


SHOULD BF CALLED 
D 
A 
B 
C 
E 


DEFINING PARAMETER GROUP N0 o t 
1 / 4 REP — ROW 1 

********************************* 
*1 *7 TEMP PRESS 

0* G MEAN 

TEMP 

*2 PRESS 

*3 TIME 

*1 43 TEMP TIME 

*4 VEL 

♦ 5 ANGLE 




4 4**4 

0.168000 

C 

C. 800000 
C. 400000 
C. 800000 
C.800000 
C. 120000 
C. 100000 


* * * * * 


DEFINING PARAMETER GROUP NO. 
1/? REP— ROWS 1,2 


4 4 4 4 4 

* 444444444 

4 4 4 4 4 4 

* * * 

* * * 

* * * 

444444444444444444444 

*i *2 


TEMP 

PRESS 



C. 280000 

0 + 


G MEAN 



C 

42 *3 

*4 


PRESS 

TIME 

VEL 

C 



TEMP 




C.800000 

*2 



PRESS 



C.800000 

43 




TIME 


C.800000 

♦ 1 43 


TEMP 


TIME 


C.8C0000 

*2 *3 



PRESS 

TIME 


C. 480000 


*5 





ANGLE C. 200000 


*4 




VEL 

1.000000 

*1 

*4 

TEMP 



VEL 

C. 500000 

*2 

*4 


PRESS 


VEL 

0 



*1 

*3 *4 *5 


*1 *2 

*4 

TEMP 

PRESS 


VEL 

C 




*3 *4 *5 


*3 

*4 



TIME 

VEL 

C. 500000 

*1 *3 

*4 

TEMP 


TIME 

VEL 

C.400000 

*1 *2 *3 

*4 

TEMP 

PRESS 

TIME 

VEL 

c 





*4 *5 


DEFINING PARAMETFR GROUP NO* 

3 





FULL— ALL 

ROVS 






4 4*44 

44*4444444 

4 4 4 4 4 4 

* * ★ 

* * * 

* * * 

444444444444444444444 

*1 *2 


TE MP 

PRESS 



0.400000 

0* 


G MEAN 



c 

*2 *3 

*4 


PRESS 

TIME 

VEL 

c 

*1 *2 ** 

*5 

TE MP 

PRESS 

TIME 


ANGLE C 

+ 1 

*4 *5 

TE MP 



VEL 

ANGLE C 

*1 


TEMP 




C.800000 

*? 



PRESS 



C.800000 

*3 




TIME 


C.800000 

*1 *3 


TEMP 


TIME 


C.800000 

42 *3 



PRESS 

TIME 


C.800000 

*1 42 *3 


TEMP 

PRESS 

TIME 


C.800000 


*4 




VEL 

1.000000 

♦1 

*4 

TEMP 



VEL 

C. 500000 


4 4 4 ** 4 4 


4 4 4 4 4 4 


32 



*1 


*2 *4 

*? *4 

*3 *4 

*1 *3 *4 

*1 *? *3 *4 

*5 

*\ *5 

*2 *5 

*1 *2 *«; 

♦3 *5 

*1 *3 *5 

*2 *3 *3 

*4 *5 
*2 *4 *5 

*1 *2 *4 *5 

*3 *4 *5 

*1 *3 *4 *5 

*2 *3 *4 *5 

*1 *2 *3 *4 *5 



PRESS 


VEL 


0 

TEMP 

PRESS 


VEL 


0 



TIME 

VEL 


C. 500000 

TEMP 


TIME 

VEL 


C. 400000 

TEMP 

PRESS 

TIME 

VEL 

ANGLE 

C 

1.000000 

TEMP 




ANGLE 

C. 400000 


PRESS 



ANGLE 

0 

TE MP 

PRESS 



ANGLE 

0 



TIME 


ANGLE 

C. 300000 

TEMP 


TIME 


ANGLE 

c 


PRESS 

TIME 


ANGLE 

c 




VEL 

ANGLE 

0 


PRESS 


VEL 

ANGLE 

c 

TEMP 

PRESS 


VEL 

ANGLE 

0 



TIME 

VEL 

ANGLE 

0 

TEMP 


TIME 

VEL 

ANGLE 

c 


PRESS 

TIME 

VEL 

ANGLE 

0 

TEMP 

PRESS 

TIME 

VEL 

ANGLE 

c 


FOR THE ABOVE PERMUTATION THE EXPECTED UTILITY IS 0.37074 


THF EXPECTFC UTILITIES AT THE 

DEFINING PARAMETER GROUP 1 

DEFINING PARAMETER GROUP 2 

DEFINING PARAMETER GROUP 3 


STQPPI NG POINTS ARE. • 
0.39850 
0. 41000 
0. 29062 


THIS MATCHING 

MAXIMIZES THE 

EXPECTED VALUE AT THE 

VAR I ABl F 

SHOULD BE 

CALLED 

1 TFMP 


C 

? PRESS 


B 

T IMF 


D 

4 VEl 


A 

5 ANGLE 


E 


2 STOPPING POINT 1/2 REP— ROWS 1,2 


DEFINING PARAMETER GROUP NO. I 
1 /4 ppp — pnw I 


♦ * 

******** 

************** 

******************** 

******* 

♦ I 

♦ 2 *3 

TEMP PRESS 

TIME 

0.200000 

0* 


G MEAN 


C 

*i 


TEMP 


C. 560000 


*2 

PRESS 


C. 400000 


*4 


VEL 

C. 200000 


*■* 


TIME 

C. 480000 


* 5 


ANGLE 

C. 200000 


*2 *3 

PRESS 

TIME 

C. 480000 


DEFINING PARAMETER GROUP NO. 
1/3 p F p — RO VS 1,2 


***** 

* *********** 

* * * * 

* * * 

* * * 

************************ 

*3 

*4 



TIME 

VEL C. 250000 

0* 


G MEAN 


C 

*1 *? 

*4 

TEMP 

PRESS 


VEL 0 

*1 


TEMP 



C. 800000 

*2 



PRESS 


C. 800000 

*1 *2 


TEMP 

PRESS 


C. 800000 

★ 3 




TIME 

C. 800000 

*1 *3 


TEMP 


TIME 

C. 600000 

*2 *3 



PRESS 

TIME 

C. 800000 

*■» *2*3 


TEMP 

PRESS 

TIME 

C. 800000 


*4 




VEL 1.000000 

*1 

*4 

TE MP 



VEL 0.500000 

*3 

*F 



TIME 

ANGLE C. 300000 

* I *3 

*4 

TEMP 


TIME 

VEL C. 400000 


*«3 




ANGLE 1.000000 

*1 

*5 

TEMP 



ANGLE C. 400000 

DEFINING PARAMETER GROUP NO. 3 





FULL — ALL 

ROVS 





***** 

************ 

* * * * 

* * * 

* * * 

************************ 

*3 

*4 



TIME 

VEL C. 250000 

0* 


G MEAN 


0 

*1 *2 

*4 

TEMP 

PRESS 


VEL C 

*2 *3 

*4 *5 


PRESS 

TIME 

VEL ANGLE C 

*1 *3 

*5 

TEMP 


TIME 

ANGLE C 

*1 


TEMP 



C. 800000 

*2 



PRESS 


C. 8 00000 

*1. *2 


TEMP 

PRESS 


C. 800000 

*3 




TIME 

C. 800000 

*1 *3 


TEMP 


TIME 

C. 800000 

*2 *3 



PRESS 

TIME 

C. 800000 

*t *2 *3 


TEMP 

PRESS 

TIME 

C. 800000 


*4 




VEL 1.000000 

*1 

*4 

TEMP 



VEL C. 500000 

*2 

*4 


PRESS 


VEL C 

*1 *3 

*4 

TEMP 


TIME 

VEL C. 400000 


********* 


33 



♦ 2 

*3 

*4 



PRESS 

TIME 

VEL 


0 

*1 

42 

★ 3 

*4 

*5 

TEMP 

PRESS 

TIME 

VEL 

ANGLE 

C 

1.000000 

*1 




*5 

TEMP 




ANGLE 

C. 400000 


*2 



4*5 


PRESS 



ANGLE 

c 

♦1 

*2 



*5 

TE HP 

PRESS 



ANGLE 

c 



*3 


45 



TIME 


ANGLE 

C. 300000 


*2 

*3 


*5 


PRESS 

TIME 


ANGLE 

c 

*1 

*2 

*3 


45 

TEMP 

PRESS 

TIME 


ANGLE 

c 




*4 

*5 




VEL 

ANGLE 

c 

*1 



*4 

45 

TEMP 



VEL 

ANGLE 

c 


*2 


♦ 4 

*5 


PRESS 


VEL 

ANGLE 

0 

*1 

*?. 


*4 

45 

TEMP 

PRESS 


VEL 

ANGLE 

c 



*3 

44 

45 



TIME 

VEL 

ANGLE 

c 

*1 


*3 

44 

*5 

TEMP 


TIME 

VEL 

ANGLE 

c 

*1 

*2 

♦ 3 

44 

45 

TE MP 

PRESS 

TIME 

VEL 

ANGLE 

c 


FOR THF ABO VF PERMUTATION THE EXPECTED UTILITY IS 0.41934 


THE EXPECTED UTILITIES AT THE STOPPING POINTS ARE.. 
DEFINING PARAMETER GROUP 1 0,31500 

OEF INING PARAMETER GROUP 2 0.59062 

DEFINING PARAMETER GROUP 3 0.29531 


THIS MATCHING MAXIMIZES THE 


EXPECTED VALUE AT THE 


3 STOPPING POINT FULL — ALL ROWS 


VAR I ABLE 

1 T FMP 

2 PPFSS 

3 TIME 

4 VFL 

F ANGLE 


SHOULD BE CALLED 
C 
D 
B 
E 
A 


DEFINING PARAMETER GROUP NO. 1 
1 /4 RFP — ROW 1 

************************************************************ 


41 *2 *3 

0 * 

*1 

*2 

*4 

*3 

*5 

*2 *3 


TEMP PRESS TIME 

G MEAN 

TEMP 

PRESS 

VEL 

TIME 

ANGLE 

PRESS TIME 


C. 200000 
C 

C. 560000 
C. 400000 
C. 200000 
C. 480000 
C. 200000 
C. 480000 


DEFINING PARAMETER GROUP NO. 2 


i n 

RFP 

ROMS 1,2 







* * 

* * * 

********* 

******* 

* * * 

* * * 

* * 

********** 

************* 


*3 

44 



TIME 

VEL 


C. 250000 

04 



G MEAN 




C 

*1 

43 

45 

TEMP 


TIME 


ANGLE 

c 

** 



TEMP 





C. 800000 


42 



PRESS 




C. 800000 

41 

4 2 


TEMP 

PRESS 




C. 800000 


43 




TIME 



C. 800000 

41 

43 


TEMP 


TIME 



C. 800000 


*2 43 



PRESS 

TIME 



C. 800000 

*1 

*2 *3 


TEMP 

PRESS 

TIME 



C.8000C0 



44 




VEL 


1.000000 

*1 


44 

TEMP 



VEL 


C. 500000 


43 

45 



TIME 


ANGLE 

C. 300000 

*1 

43 

44 

TE MP 


TIME 

VEL 


C.400000 



*5 





ANGLE 

1.000000 

41 


45 

TEMP 




ANGLE 

C.400000 


DEFINING PARAMETER GROUP NO. 3 


FULt 

— 1 

ALL 

ROMS 





* * 

* > 

*> * 

* ******** 

******* 

* * * 

* * * 

* * 


42 


45 


PRESS 



0* 




G MEAN 



41 


*3 

45 

TEMP 


TIME 



42 

43 

44 45 


PRESS 

TIME 

VEL 

*1 

*2 


44 

TEMP 

PRESS 


VEL 





TEMP 





*2 




PRESS 



*1 

42 



TEMP 

PRESS 





43 




TIME 


*1 


*3 


TEMP 


TIME 



*2 

43 



PRESS 

TIME 


41 

42 

43 


TEMP 

PRESS 

TIME 





♦4 




VEL 

41 



44 

TEMP 



VEL 


42 


44 


PRESS 


VEL 



43 

44 



TIME 

VEL 


****************************** 

ANGLE C 

C 

ANGL E C 

ANGLE C 

C 

C. 800000 
C. 800000 
C. 800000 
C. 800000 
C. 800000 
C. 800000 
C. 800000 
1.000000 
0.500000 

c 

C. 500000 
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I 


*1 


*3 

♦ 4 


TEMP 


TIME 

VEL 


C. 400000 


*2 

*3 

*4 



PRESS 

TIME 

VEL 


C 

*1 

*2 

*3 

+ 4 

*5 

TEMP 

PRESS 

TIME 

VEL 

ANGLE 

c 

1.000000 

*1 




*5 

TEMP 




ANGLE 

C. 400000 

*1 

*2 



*5 

TEMP 

PRESS 



ANGLE 

0 



*3 


*5 



TIME 


ANGLE 

C. 300000 


*2 

*3 


*5 


PRESS 

TIME 


ANGLE 

c 

*1 

*2 

*3 


*5 

TEMP 

PRESS 

TIME 


ANGLE 

c 




*4 

*5 




VEL 

ANGLE 

c 

*1 



*4 

*5 

TEMP 



VEL 

ANGLE 

c 


*2 


*4 

*5 


PRESS 


VEL 

ANGLE 

0 

*1 

*2 


*4 

*5 

TEMP 

PRESS 


VEL 

ANGLE 

c 



*3 

♦ 4 

*5 



TIME 

VEL 

ANGLE 

c 

*1 


*•* 

*4 

* 5 

TEMP 


TIME 

VEL 

ANGLE 

c 

*1 

*7 

*3 

*4 

*5 

TEMP 

PRESS 

TIME 

VEL 

ANGLE 

c 


FDR THF ARHVE PERMUTATION THE EXPECTED UTILITY IS 0.42169 
THE EXPECTEC UTILITIES AT THE STOPPING POINTS ARE.. 


DEF INING 

PAPAMFTER 

GROUP 

1 

0. 31500 

D C F INTNG 

PARAMETER 

GROUP 

2 

0© 59062 

defining 

PARAMETER 

GROUP 

3 

0.30312 


SUMMARY OUTPUT T ABL E 




** BAYES *** 

**** x ***** 

*+** 2 ***** 

**** 3 

C 

0 

C 

C 

c 

A 

B 

0 

e 

B 

0 

B 

E 

C 

A 

E 

A 

E 

E 

A 


EXPFCTED UTILITY 
OVFR STOPPING PTS 


0 42169 


0. 37074 0.41934 0.42169 


EXPFCTED UTILITY 
AT EACH STOPPING P 
* 0 © 3 1 *00 

? 0^9C62 

3 0.30312 

CURRENT EXECUTION 


0© 39850 
0. 41000 
0. 29062 
TIME 


0.31^00 
0. 59062 
0.29531 

0© 07 


0.31500 

0.59062 

0.30312 
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APPENDIX B 


FORTRAN LISTING 


*!0FTC BLOK 

BLOCK DATS 

COMMON/BDATA/ POWERSUl), I ALPHA <9 ) , ALPHAUO), 

X IUNIN, IUNOUT, MASK, NEG 

C* ************************** 

IMTFGFR LOWERS 

DATA (POWERS (11,1=1,11 ) /0 , 1 ,2 ,4, 8, 16 , 32 ,64 ,12 8 , 256 , 0/ 

DATA ( ALPHAU) ,1=1,10)/ 1 HA , 1HB ,1HC , 1H0, 1HE, 1HF ,1 HG, 1HH, 1 HJ , 1H / 

DATA TUNTN/5/, IUNOUT/6 /,MA SK /01 / »NEG/G4000000 00000/ 

END 


* I PFTC NAMER DFBUG 

COMMON/BOA TA / POWERS(U), IALPHA(9), ALPHA! 10), 

X TUN IN , IUNOUT, MASK, NEG 

COMMON /Bl./ NFAC,NCLASS,NN,NDPG,PBAYES,IRUN,PBAYEX, IUT ILF, UTSWCH 
COMMON /BP/ 0PG(128,32) , BLOCK (128 ,32 ) ,I0PG(32> , NBLOCM 3 2 ) , I P( 9 ) , 

X ID< 9),NSURl (9) ,LPFRM< 9) ,1 I (10) ,PSTOP(32) , WT (3 2 » , PSU MX < 32 ) , PSUMC 32) 
COMMON /R F ST/PR OR ( SI 2) , Tl(512), PBL0CK(128 ,32 J , IS AVER ( 9, 3 3 ) , 

X UTIL(fl2>, UL I ST ( 1 28 ) , IULIST(128), ID0CG(4,32) 

INTEGER POWERS, TJ. , DPG, BLOCK 
LOGICAL UTSWCH, LPERM 


PEAL TDENT 1 C 

EQUIVALENCE (X,!X) 11 

DIMENSION DUMPl(l) , 0 UMP 2 ( 1 ) 12 

EQUIVALENCE (DUMP! ( 1 ) ,NFAC ) , (DUMP 2 I 1 ) , IDPG(l) ) 13 

DIMENSION IT 1 ( 8 ) , XT 3 ( 9 ) , I DENT ( 1 A) 14 

DATA BLANK / 6 H /, ENDCR.D/ 6 HENDALL/ 15 

COMMON /L INX/ XN 0 UT( 5 ) ,H 0 L 0 UT( 9 ) ,FACU 4 , 9 > 16 

C 17 

C* ********************* ****************************** **** ******* ******** 18 

LQ 1 = 9 19 

CALL T IMF 1 (TSTART) 2 C 3 

TMAX= C.O 21 

LQ 2 = 238 22 

23 

* ** ******* ******** ***************** *************** ************* ******** 2 A 

25 

NAMER SECTION 1 26 

27 

1.0 R FA0( IUN IN * 5000 ) IDENT 28 5 

IF( IDENTd. J.EQ.ENDCRD) STOP 29 

WR ITF( IUNOUT, 5005) TDENT 30 

31 

****** ******************** ************* ************* ********* ********** 32 

33 

NAM FR SECTION 2 34 

35 8 

DO 12 J = !» 5 36 

12 XNOUT ( J ) =8LANK 37 

no 14 J— 1 , 9 38 

I ALPHA ( J 1= J 35 


36 


vOODJfMn^UMM ^00«JOU)>(v)MH 



o d n "i n o -» o “> o .n n r> o -i r» n o r> o n 


1 A HOinUT( J 1=RLANK 

READ< IUNIN, 50101 TMAXX 
TMAX= TM AX + TMAXX 
WR ITF( IUNOUT, 50701 TMAXX 
CALL T IMFKTPR TNT) 

TPR IN T = ( TPRINT-TSTART) /3600.0 
WP TTF( IUNOUT, 50251 TPR I NT 

*********** **************** ******************************************** 

NAM FR SECTION 3 

REAO( IliNIN,50451 TTYPRN 
GO TO <30,154,?01,ITYPRN 
20 RFADt IUNTN, 50701 ( XT3 ( K ) , K=1 ,9 1 

* ********* * ************************** ********************************** 

NAMER SECTION A 

30 REAO( I UNIN, 5045) NF AC 

IF( (NFAC.l T.l) .OR. (NFAC.GT.91) GO TO 8020 

WP TTF( IUNOUT, 50501 NF AC 

I Ff ITYPRN.NF.31 GO TO 38 

on 35 K = 1 fNFAC 

DO 3* L = 1 , N r A C 

LL=t 

I<MXT3(K laFQoALPHA(L) 1 GO TO 3A 
3 3 CONTINUE 
go Tn PO^ 0 

° A IGAVEP (K,1 ) = L L 
7* CONTINUE 
7p CON tin U c 

*********************** ********** **** *** ********** ********************* 
N AM p P SFCTION 5 
on AO J = i f NFAC 

RFAD( I UN IN, 50 00) < FAC ( I , J 1 , 1 =1 ,141 
WP ITFC TUNOUT, 50 551 J, (F AC (I , J ) , I =1 ,1 A 1 
^0 CONTINUE 

*********************************************************************** 
N AM FR SFCTION 5 
TTf 1 1 = ^ 

l P FPM ( i )=, TRlJFo 

NSUB! ( 1 1=NFAC 

REAO( TUNINfF045) NCLASS 

TF( (NCLASS.LTo 1 1o ORo (NCLASSn GT, 9)1 GO TO 8030 
IFfNCL ASS-11 60,60,50 
50 p FAOf T IJN IN , 5045 ) ( NSUB T (I 1,1=1, NCLASS) 

WR T Tp f IUNOUT, 7000 1 NCLASS, (NSUB I (I) ,1=1 , NCLASS) 

OH 55 J=1,NCLA$S 
L P FRM ( J 1= Q TRUF n 
I T ( J + 1 1 = I T ( J 1 +NSUB I ( J 1 
5 5 CONTINIF 
A 0 NN= 3*+NFAC 
nn 65 J=1 T NN 
TM J1=c 

PPOR( J ) = “* 0 

UTIMJ 1 = 0 0 
f> c CONTINUE 

******* ********************************** ********* ********************* 


AC 


A 1 

24 

A2 


A3 

25 

AA 

26 

A 5 


A6 


A 7 


A 8 


45 


5C 


51 

28 

52 

29 

53 


5 A 


55 


56 


57 


56 


55 

32 

6C 

35 

61 


62 

A3 

63 


6 A 


65 


66 


67 


66 


65 


7 C 


71 


72 


73 


7 A 


75 


76 


77 


78 


75 

65 

ec 

7A 

81 


82 


83 


8 A 


85 


86 


87 


88 


85 


90 

e2 

91 


92 


53 

88 

9 A 

55 

95 


56 


97 


98 


55 

113 

IOC 


101 


102 


1 C 3 


1C A 


105 


106 



37 


r 




107 


c 


NAMFR SFC.T ION 7 


ioe 


r 




105 




RFADI I l!N IN, 5048) NPIN.IUTILF ,UC0EF 


lie 

123 



WRITE ( IUN0UT»6045) N° I N» I UTI LF 


111 

126 



utswch= 0 falsf,. 


112 




IF( I IUTH f.LT»1I „OR, (IUTILF.GT.5) ) GOTO 

8035 

113 




TFUNP IM,LF, OI»OR„ (NR INoGT.NN) ) GO TO 8032 


114 




!PUUrnFF,lTo0.0» o 0R> (UCOFFoGT.1.0)) GO TO 8033 

1) 5 




TPf IUTIL^GF.R 1 UT55wrH= a TRUFo 


116 




TP( IUT Tl^EO, 5) WR !TE ( T tJNOUT ,7010 1 UCOEF 


117 

138 



nn PO J=1,N°TN 


118 




R EAH( T UN IN , 50 60 1 C in ( I ) , T =1 , 9) , P, UT 


115 

142 



TPf (P oCT JcOlflOR, (Pet 1*0, 0) 1 GO TO 8038 


12C 




T =0 


121 




T T = 0 


122 




KK =0 


123 




on 67 K = i,9 


124 




K K X - If-K 


125 




K 1= ITI(KKX) + 1 


126 




T = T + Pnw8RS(KT) 


127 




TFCKT-J) 67,67,66 


128 



66 

TT = IT 4 ITM KKX) *1 0**KK 


125 

162 



KK =KK + 1 


13G 



6 7 

TON T IN UE 


131 




1= 1 + 1 


132 



8 * 

OH to (P?, R4, 86,99 ,87) ,IUTTLF 


133 



8 ? 

UTIL( I 1 = 1,0 


134 




GO TO P9 


135 



8* 

UTIL( T 1=P 


136 




GO TO F9 


137 



RA 

UTIt f I 1=UT 


138 




GO TO 39 


135 



R 7 

OT I L ( T )= UCnFF*UT + ( 1 0 0- UC OF F ) * P 


14C 




On to pq 


141 



R 8 

«JTTL< I 1=UT*P 


142 



R 9 

continif 


143 




WPTTF( I UNO UT , 6060 ) I T , P , NT I L < I 1 


144 

182 



PP 0 B( T 1 = 1 « 0 - P 


145 



QO 

CON T I N UF 


146 


C 




147 


c *** ******* * *********** ***** *********************** ********************* 

148 


c 




149 


r 


NAM5R SECTION 8 


X5C 


r 




151 




PEAnC I IN IN , 504 5 ) NOPG 


152 

187 



IF( (NDPG,l_T>1 )oOR, (N0PGoGTo32)) GO TO 8140 


153 




WP ITE( TUNOUT, 60651 ND PG 


154 

191 



on 180 1=1, NOPG 


155 




PE ADI I UN IN, 50001 ( I0OCG(K,I) ,K=l ,4) 


156 

154 



WR T TF ( IUNPUT, 6 0641 (IODCG (K,I 1 ,K=1,4) 


157 

199 



PFADI I UN IN, 3 06 31 NGFN,PSTOP(I) ,WT(I) 


158 

204 



IP( (NG c NoLTo 01 oQP t (NGFNoGToMTN0(NFAC,71) 1 GO 

TO 9150 

159 




IF( (PSTPPI IlolToO.OIoORo < P STOP <I>.GT. 1.0) ) GO TO 8160 

16 C 




WP TTF( TijNOUT, 6C66) I,NGEN,PSTOP (11 , WT ( I ) 


161 

215 



I OP G( I 1= 2**NGFN-1 


162 


r 




163 


c* 

* 5{ 

*********************** 

* * 

164 


r 




165 


r 


NAM r R SFCTION 3A 


166 


r 




167 

218 



IF(NGFN) 106,106,91 


168 



qi 

no 103 j=i ,ngfn 


169 




REAO( I UN IN , 50 7 0 ) ( XT3 ( K) , K=J. ,91 


17C 

223 



WP TTE( IUNni)T,6070) (XT->(K) ,K=1,9) 


171 

228 



K T = 0 


172 




90 JOO K = 1 , 9 


173 



38 





nn q4 L=1 ,T 0 

174 




LL =L 

17 5 




T f ( XT 3 ?K ) 0 fq ALPHA ( LH GO TO 98 

176 



9* 

CONTINUE 

177 




go T n po^o 

178 



o p 

KKT=LL+1 

178 




KI = K I + PnwFRS(KKI) 

18C 



1 00 

CONTTMUP 

181 




nor y ( j, i )= KI 

182 



i 0 *? 

CONTINUE 

183 



i Of 

ioi= iopnrn 

184 




IF? TOI-1 ) 1.21 ,1?1 ,108 

185 



TOP 

KP = 4 

186 




M1= NGFN 

187 




LPT'RstNJGFN + I 

lee 




LL PM = 0 

188 




DO 120 J=? ,NGFN 

18C 




KK = J— t 

191 




iF(npr,(j f n,F 0 o DPGn »in go to 806 o 

192 




on 110 K=1 ,KK 

193 




DPG(LPTR t T)=IEXOR(DPG( J,I) ,0PG(K,I)) 

184 

275 



IF? DPGd PTR, I ) o F0 OPGUvin GO TO 8060 

185 



V ! 0 

lptr =lptr + i 

186 




T F? L L FA© PO oO ) Gn TO 118 

187 




DO 1 1 8 K = 1 ,LL FN 

198 




N=N GFN 4K 

188 




PPG? LPTR , I )=IFXOR (DPG( J, I ) ,DPG(N,I ) ) 

2CC 

282 



IF(DPG(LPTR f M,FQ„0PG(1 ,T)> GO TO 3060 

201 



V c 

I ptr=lptr+i 

202 



1 1 p 

LLFN=?*LLFN+KK 

203 



* ->o 

fONTTNlIF 

204 


r 



205 


r* 

: * j 

************************* 

206 


c 



207 


r 


NAMFR SFCTION 8B 

206 


r 



208 



1 ? 1 

RFAD? HJNIN,5045) NBLOCKCT) 

210 

3C7 



NR = N BLOCK? I ) 

211 




TFf (NP ,1 T,0),PPo (NRoGT«NN) ) GO TO 81.70 

212 




WR !TF( IUNOUT, 6080) NR 

213 

313 



IF(NB) tcq f'SOtl?.? 

214 



1 9 ? 

on 140 J=i ,NR 

215 




RFAD? ILNTN,5070) (XT3(K) ,K=1 , 9 ) , PR LOG K ( J , I ) 

216 

317 



IF? (PBLOf.KU, IJeLToO, 0)o ORo ( P8L0CK ( J f I) • GT* 1, 0 ) ) GO TO 818 0 

217 




WR ITF? HJNOUT » 60 7 0 ) ( XT3 ( K ) ,K=1 ,9) ,PBLOCK( J ,1) 

218 

327 



c 

ii 

* 

218 




no 130 K =1 t 9 

22C 




no 194 1 =1,1 0 

221 




LL =1 

222 




IF? XT3 (K ) o F0n ALPHA? L) ) GO TO 128 

223 



1 74 

CONTINUE 

224 




Gn TO fO^O 

225 



’ 9 R 

KK T=LL4i 

226 




KI = K I + PnWFR S( KKT ) 

227 



** 30 

rPNTINUF 

228 




BLOCK? J, T) = KI 

228 



1 4 0 

CONTTMUF 

230 . 



5 r f) 

contincf 

231 




IF? TTY D PN^ FQo 3 ) CALL 0NCF($1Q) 

232 


r 



233 


0* ****** *** ******* ********** ** ***************************** ***** ******** 

234 


r 



235 


r 


namtr ^fction 9 

236 


r 



237 

359 



dbayFS=C«0 

238 




nn 182 I=i,NDPG 

238 




pgum? I ) = 0 0 

24 C 



39 



16? CnNTlNUF 



241 


TPIJN = 7Q 



242 


CALL T TMF1 CTPR TNT) 



243 

370 

TPRINT = (TPRINT-TSTART) /3600a 0 



244 


WP IT£f IUNOUT, 5025) TPRTNT 



245 

372 

on to is* 



246 


r 



247 


C*** *************** ********* *********************** ********************* 

248 


r 



249 


C NAM ER SFC T ION 10 



25C 


r. 



251 


154 CONTINUE 



252 


WR ITF ( IUNOUT, 5040) 



253 

374 

CALL BCREADt DUMP1 ( L ) » 0UMP1 ( LOT ) ) 



254 

376 

CALL 9CRFA0(DUMP?<1 ) , DUMP ? (L02 ) ) 



255 

379 

CALL BCREADI IALPHAU) , I ALPH Af NFAC) 1 



256 

382 

CALL BCREA0(FACU,1),FAC (14, NFAC)) 



257 

385 

CALI 8CRFAD ( PR08 ( 1 ) , PROB ( NN ) ) 



258 

388 

CALL .8 CR EAO( UT IL ( 1 ),UTTL(NN) ) 



255 

351 

CALI 8CREAD <1 SA VFP (1 , 1 ) , I SA VEP < NF AC , 1 ) ) 



26C 

354 

CALL BCR EACH I DDCG f X ) ,IDDCG (4, NDPG) ) 



261 

357 

00 156 I = \ * NDPG 



262 


T 0 1 = TDPG ( I ) 



263 


N8=NBLCCK( I) 



264 


IF( ID I .NFo 0) CALL 8CRE AD (OPG < 1 , I ) ,OPG(IOI,I)) 



265 

409 

CALI 8CREA0 ( I SA VEP (1 , 1 +1 ) , I SAVE P( NFAC , 1 + 1 ) ) 



266 

413 

IF ( NBa EQ„ 0 ) GO TO 156 



267 


CALL BCREA0(PBL0CK(1,I ),PBLOCK(NB ,t)( 



268 

420 

CALL BCREAD(BLOCK(L ,1) ,BLOCK(NB ,11) 



269 

424 

156 CONTINUE 



27C 


c 



271 


C ******************* ******* ************************************* ******** 

272 


1 *5 R CONTINUE 



273 


r 



274 


r ****************** ********************************************** ******* 

275 


r 



276 


C NAM FR SECTION U 



277 


C EVALUATE THE CURRENT ORDERING 



278 


C 



275 


DO 100C 1 = 1 * NDPG 



28C 


r 



281 


Q* **********< s ***** 



282 


c 



283 


C INITIALIZATIONS 



284 


r 



285 


I RUN= I RUN + l 



286 


PSIJMY ( I ) =0 o 0 



287 


NB=NBLCCK( I ) 



288 


101 = inPG(T) 



289 


1011= I rr*i 



29C 


C 



291 


C*-.* ,*_*_*-*-*_*_*-*-*- *-*_*-.*_*_*-*_* _*-.*_*_*_*_*_ 

*-*-* 

* 

i 

* 

i 

* 

i 

* 

l 

* 

l 

* 

l 

* 

i 

252 


r 



293 


C NAM ER SECTION 11 A 



254 


c 



295 


C CHECK FOR BLOCK CONFOUNDING^ IF THERE ARE NO 

BLOCK 

PARAMETERS 

256 


C GO to NAMER SECTION 1TB 



257 


C 



258 


165 TPfNR) 228,228,166 



255 


ia* 00 2?6 K=l,NB 



3CC 


K 1= BLOCK ( K , I ) 



301 


KP = KT+1 



3C2 


T1(KI1 )= TRUM 



303 


C 



3C4 

■joc 







C 



3C6 


C IS this A FULL OR FRACTIONAL FACTORIAL 



3C7 



40 



o o n i n -> n r> -> ^ n -> n o o r> 


r ace 

I FI TO I ) 7035, 167,1 71 309 

’67 PSUMX ( I ) = PSUMXm+UTXL(KTlt*(l, 0-PBLOCK(K,I) » 310 

on to ??6 311 

312 

IF BLOCK PAR A METER HAS PRIOR PROB=1.oO THFRF CAN RE NO UTILITY 313 

FOR THIS ESTIMATOR,, JUST TAG ALIASES,, 314 

315 

VM TFtPBLOCKIK, Il-loO) 1.76,173,173 316 

17-» no 174 TK= ’ , IOI 317 

JJ = I EXOP. <KI ,DPG< IK, I n+1 318 461 

T 1 1 J J > = IR UN 319 

174 FONT IN L'F 320 

GO to ?P6 321 

322 

323 

324 

IN ITT AL T7E BEFORE FINDING OPTIMAL MATCHING OF PARAMETER TO ESTIMAT 325 

326 

PS= PRCRIK II ) 327 

PR= PS 328 

TUL 1ST ( 1 )= K 1 1 329 

UL I ST ( 1 )= U T I L ( K II ) 330 

T ST AP = C 331 

TF(o<;) 178,1 78,179 332 

’78 T STAR = 1 333 

Ut I ST ( 1 ) = -IJTILIKI1 ) 334 

PR=’ 0 335 

17<5 TMA X = KT1 336 

337 

338 

339 

FOR UTILITY FUNCTIONS 1 AND 2, PARAMETER WITH MAXIMUM PROBABILITY 340 

HAS MAXIMUM UTILITY, 341 

342 

I F( UTSWC H ) GO TO 200 343 

00 ’87 KK =’ , I D I 344 

JJ = TFXOR <KI ,OPG(KK, I l)H 345 484 

T1(JJ)= IR UN 346 

IF(PPne<JJJ) 2035, 1 84 ,18'* 347 

1.97 PR = PR*PRnB{JJ) 348 

1 F( PROP ( J J )-P s ) 1 86,1 86,187 349 

1PA I ST AR = ISTAR+’ 350 

* 86 P S= PROB( JJ) 351 

IMAX= JJ 352 

i 8 7 CONTINLF 353 

354 

COMPIJTF UTILITY 355 

356 

TF( ISTAR-1 ) 190,1.92,276 357 

’90 PR= PR /PROB( IMAX) 358 

l.o? PSUMXl I)=PSUMX( I M-PR*UTIL ( IMAX)* t]. 0 O-PBLOCMK, I) ) 359 

GO TO ??6 360 

361 

362 

363 

FOR ARBITRARY UTILITY FUNCTIONS COMPUTE UTILITY FOR EACH MATCHING 364 

365 

700 DO 706 KK = 2» ID I’ 366 

J J= TFXOR (K I,D»G(KK-1., I ) ) + l 367 511 

T1(JJ)= IR UN 368 

IUL I ST ( KK ) = JJ 369 

IF(PROBIJJ)) 704,204,702 370 

?0? PR= PR * PPOB(JJ) 371 

UL TSTIKK )= IJTIL(JJ) 372 

GO to ?06 373 


41 


(J L 


II 


90^ TSTAP = I STAR + 1 374 

UL TSTIKK ) = - UT I L ( J J ) 375 

906 CONTINUE 376 

377 

FIND MATCHING THAT MAXIMIZES UTILITY 378 

379 

IF( ISTAR-i ) 208,214,296 38C 

20 8 on ?1C K K = 1 , inn 381 

JJ= IIJLIST(KK) 382 

UL T ST ( KK ) = UL T ST< KK )* PR /PR08 < JJ) 383 

9^0 CONTINUE 384 

IJMAX-O oO 385 

On ? KK=] ,1011 386 

IF(ULI C T I KK ) - UMAX ) 212,212,211 387 

91 1 UM AX = ULISTfKK ) 388 

? CONTTNUF 389 

P SUMX ( T ) = PSUMXm+UMAX*( 1.0 - PBLOCK ( K, I ) ) 39C 

GO TO ?26 391 

C 392 

C 393 

91 a On 91 6 KK = I, 1011 394 

K S = KK 395 

TST = AND ( UL I ST ( KK ) , NEG I 396 

T ST = OR (TST, MASK) 397 

!F( TST ) 9 t 2035, 216 398 

21 6 CONTINUE 399 

2 1 R P S'JMX ( I ) - PSUMX(I) - ULI ST ( K S ) *PR* C lo 0-P8L0CK ( K , I ) ) 40C 

?9A CONTINUE 401 

C 4C2 

C* ***************** 403 

r 404 

2?p IF(IOI) 229, 229,245 405 

C 40 6 

r *_*_*_*_*_*_*_*_*_*_*_*_*_*_*_*_*_*_*-*_*-*..*_ *_*_*-*_*_*_ *_*_*_*_*_* 407 

C 408 

C NAMER SECTION UB 409 

C BLOCK PARAMETERS HAVE NOW BFEN ACCOUNTED FOR. 41C 

C CONTINUE COMPUTING UTILITY FOR REMAINING PARAMETERS,, 411 

C COMPUTE UTILITY EOR FULL FACTORIAL 412 

C 413 

299 IF( HIT ILFoNFol) GO TO 93 0 414 

PSIJMX ( I ) = PSUMXm+FLOAT(NN-NBLOCK(in 415 

GO TO 10CO 416 

230 nn 934 K=!,NN 417 

IF(THK)-TRUN) 231,234,231 418 

23’ P S UM X ( I ) -P SUM X(I)+UTIL(K) 419 

934 CONTINUF 420 

GO to 1000 421 

C 422 

C 423 

C COMPUTF UTILITY FOR FRACTIONS 424 

C 425 

C FIND NEXT UNTAGGED PARAMETER 426 

r 427 

245 DO 7C0 K=1,NN 428 

IF(TM Kl-TRUN) 246,700,246 429 

9 *6 PS= PROB(K) 43C 

KM 1 = K-1 431 

PP=PS 432 

TSTAP=C 433 

IUL TST(1 )= K 434 

UL I ST ( } ) = UTIL(K) 435 

IF(PS) 248 , ?4ft, 950 436 

I ST AP = 1 437 

UL T ST ( 1 )= -UTTL(K) 43e 

PR=^0 439 


42 


I 



-> r> n r> r> -> n o n o "> r> o n n r> o n 


750 T M AX=K 44C 

IF! UTS WCH ) GO TO 430 441 

r 442 

443 

444 

for UTILITY FUNCTIONS 1 AND ?, PARAMETER WITH MAXIMUM PROBABILITY 445 

HAS MAXIMUM UTILITY,, 446 

447 

DO A10 KK =1 » I D I 448 

JJ = 1 EXORIKM1 ,OPG(KK , I ))+l 449 609 

TUJJ) = TPUN 450 

IFIPROPIJJH 380,390,380 451 

380 PR = PR *PROB( J J ) 452 

IF( PRORI JJ )-PS I 400,410,410 453 

390 T STAR = T STAP+I 454 

400 PS=PROB( JJ ) 455 

TMA X = J J 456 

4’ 0 CONTINUE 457 

GO to =05 458 

455 

4 6 C 

461 

FOR ARBITRARY UTILITY FUNCTIONS COMPUTE UTILITY FOR EACH MATCHING 462 

463 

430 DO 460 KK=2,IDIl 464 

JJ= IEX0RIKM1 ,DPG(KK-1 ,1 ) ) +1 465 628 

TMJJ) = IP. UN 466 

ItlL 1 ST ( KK ) = JJ 467 

IF ( PROBI J J ) I 450,450,440 468 

440 PR= PR*PPOB(JJ) 469 

IJL T ST I KK 1= UTIL ( J J I 470 

GO to 460 471 

450 I ST AR = I STAR +1 472 

ULTST(KK!= -UTILIJJI 473 

440 CONTINUE 474 

GO TO 600 475 

476 

477 

478 

INCREMFNT UTILITY OF THIS DCG BY UTILITY CF ESTIMATOR 475 

480 

505 IF(ISTAR-I) 520,530,700 481 

5?C PR= PR /PROBI IMAX) 482 

530 PSUMXI I)=PSUMXm+PR*UTIL(IMAX» 483 

GO td 700 484 

485 

486 

487 

600 IF(ISTAR-l) 620,660,700 488 

4? 0 DO 630 KK=!,I0I1 485 

JJ= I UL I ST IKK ) 490 

IIL I ST ( KK I = Ul. I ST I KK I * PR/PROBIJJ) 491 

67 0 CONTINUE 492 

|JMAX = 0 ,0 493 

DO 640 KK = 1.» I D II 494 

I F ( UL I STIKKI-UMAX1 640,640,635 495 

UMAX= UL 1ST! KK I 496 

6 40 CONTINUE 497 

PSUMXI I»= PSUMXI I l+UMAX 498 

GO tq 70C 455 

660 DO 670 KK=1,IDIT 50C 

K5= KK 501 

tst=AND( UL 1ST! KK) , NFG I 502 

TST=OR (TST.MASK) 503 

TFITST) 680,2035,670 504 

670 CONTINUE 505 

690 PSUMXI I)= PSUMXI!) - ULIST(KS)*PR 506 


43 


o r> n r> n o n n -i r» n o r> 


^00 continuf 
1 000 continuf 

C 

f ************************************* ** ******* ****** ******************* 

c 

C NAM FR SECTION 12 

r 

P B A YF X = 0. 0 

on 10?C 1=1 1 NDPG 

° SUMY { I)= PSUM X { I ) *WT ( I ) 

P B A YF X = PBAYFX + P SUM X (I) * PSTOP (I) 

10?0 COMTINUP 

IF( PBAYFX-PRAYES) 1060,1060,1030 
1 030 P BA YF S =P BA YFX 

On 1050 T=1,NFAC 
I SA VEP ( T , 1 ) = I ALPHA ( I ) 

1 050 CONTINUF 

10^0 on It cc 1 = 1, NDPG 

I F( PSUMX ( T )— P SUM ( I ) 1 1.100,1100,7 070 
1 070 PSUMf I ) =P SUM X ( I ) 
nn 1080 J=1,NFAC 
TSAVFP (J, I+U= I ALPHA ( J ) 

1 0*0 CONT IN UP 
1 i 00 CONTINUE 

*********** ************************************************************ 
NAMFR section 

00 11 5 C K = ! ,NCLASS 

1 SFND = 1 1 ( K ) 

CALL PFP.MUT( IALPHAI TSENO) ,NSUBI (K ) , LPERM ( K ),ISEND) 

IF(LPPRMfK)) GO TO 1150 
GO Tn 2000 
1150 CONTINUF 

GO TO 2000 

*********************************************************************** 

NAMPR SECTION 14 

7Q00 CALL T IM F i ( TNOW) 

T = { TNOW-TSTAPTI /3^00, 

IF(T-TMAX) 7 58,2010,201 0 

*********************************************************************** 

NAMFR SECTION 15 

70" 0 CONTINUE 

W D ITF( TUNOUT, 60901 T, T RUN 

CALL BCD IJM P ( 0 U M P 1 (1),0UMP1 ( LD1 ) ,0) 

CALL BCOUM o ( HUMP? (X), DUMP? (L02) ,0) 

CALL BCDUMP ( I ALPHA (1) , I ALPH A ( NF AC) ,0) 

CALL BCOUMP(FAC (! ,1) ,FAC ( 14 , NF AC ) ,0) 

CALL BCOUMPI PPOR(I), PR08<NN>,0> 

CALL BCOUMP(UTIL (7 ) ,UTTL(NN) ,0) 

CALL RCDUMP ( TSAVEPfi tilt ISAVEPCNFACtll *0) 

CALL BCDUMPI IDOCGH ,7 ) ,IODCG (4, NDPG) ,0) 

DO 20 3 C I =1 , NDPG 
ID T= I OPGf I) 

NB=N8L OCK ( M 

TF(IOI.MPoO) CALL BC0UMP(0PGU t I) f DPG(IDI,n,0) 

CALL BCOUMP ( I SA VEP ( 1 , I +1 ) , I SAVE P( NF AC » I + l ) ,0 ) 

IF(NR o FQ o 0) GO TO 2030 

CALL BCDUMP(PBL0CK(1, I) ,PBLOCK(NB ,I),0) 

CALL BCDUM P ( BL OCK ( 1 , I ) ,BLOCK(NB ,I),0) 

20^0 CONTINUF 
2075 STOP 


507 
5C6 
505 
51C 
511 
5 J 2 

513 

514 

515 

516 

517 

518 
515 
52C 

521 

522 

523 

524 

525 

526 

527 

528 
525 
53C 

531 

532 

533 

534 

535 

536 

537 

538 742 
535 

54C 

541 

542 

543 


544 


545 


546 


547 


548 

755 

545 


55C 


551 


552 


553 


5 54 


555 


5 56 

758 

557 

760 

558 

763 

559 

766 

56C 

769 

561 

772 

562 

775 

563 

778 

564 

701 

565 


566 


567 


568 

753 

565 

797 

570 


571 

804 

572 

808 

573 


574 
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n n n n r> •**> n n n -> 


*********************************************************************** 

NAM FR SECTION 1.6 

1000 CONTINUE 

CALL T TME1 (TPR INT) 

TPRINT=(TPRINT-TSTART) /3600>0 
WP TTF( IUNOUT, 8025 > TPR T NT 
CALL REMACH 
CALL T 1MEI (TPR TN T ) 

TPR IN T = ( TPR IN T-T START) /3600. 0 
WR ITF( IUNOUT, SOPS) TPR INT 
GO Tf) 10 

*********************************************************************** 
5000 FORMAT! J.1A6»1 A?) 

5005 FORMAT (1.4H1NAMFR OUTPUT I3A6,1A2/1H ) 

50’0 FORMAT ( F6. 0) 

5020 FORMAT ( 50H PROGRAM WTLL OUMP FOR RESTART IF NOT FINISHED IN F10.0, 
X PH MINUTES. /IH ) 

5025 FORMAT ( ? 5H CURRENT EXECUTION TIME F12.2) 

5010 format (LI) 

5040 FORMATI38H THIS IS A RESTART OF A PREVIOUS CASE /I H ) 

■50A-* format (2T6»F10»9) 

5045 FORMAT (916) 

5050 FORMAT ( 1. 1.H THERE ARE I4,22H FACTORS. THEY ARE... /IH ) 

5055 fopmaT(1x,I1.0,3X,13A6,A2> 

5 06 0 FORMAT (PIT, 2FI0.0) 

5065 FORMAT! I 6, F6„ 6 , F6. 0 ) 

5070 FORMAT (SA1 ,F5 0 5) 

60^-5 FORMA T ( i HK I6,6?.H PARAMETERS WITH NON-ZERO PRIOR PROBABILITIES AND 
X UTILITIES /1X/19H UTILITY FUNCTION 17) 

6060 format ( lx, 19, 2GI4. 6) 

6065 FORMATOHR 16 , 26H DEFINING PARAMETER GROUPS /IX) 

6064 FORMAT (1HK 1 20( 1H- ) /1HK 4A6) 

6066 F0PMAT(4HK0PG I3/4H I3,11H GENERATORS / 

X ?\H PROB OF STOPPING F10.5,13H WEIGHT F10.6) 

60 T 0 FORMAT ( 5X, 9A1, FT0.5) 

6080 FORMAT (IOHKTHERE ARE 14, 1 7H BLCCK PARAMETERS /IX) 

6090 FORMAT ( 3 5H TIME EXCEEDED. DUMPING FOR RESTART /13H EXEC. TIME 
X F7 - 2 , 5H MIN. /8H IRUM = 115) 

7000 FORMAT (1IHKTHERE 4RE T1,21H CLASSES OF VARIABLES /5X.9I6) 

701 0 FORM ATI 9HK UCOEF= F1.2.9) 

*** ******* * ************************************************************ 
FRROR MF5 SAGES 

5010 WR ITF( IUNOUT, 90T0) 

9010 FORMAT ( 40H ILLEGAL CHARACTER IN SPECIFIED MATCHING ) 

GO T n 2015 

5020 WR TTE( IUNOUT, 9020) NF AC 

°0’0 FORMAT ( 33H NUMBER OF FACTORS OUT OF RANGE 16) 

GO TO 20’5 

8010 WR ITE( IUNOUT, 9010) NCLASS 

9010 FORMAT ( ?°H NUMBER CLASSES OUT OF RANGE 16) 

GO TO 2035 

8032 WR I TE( IUNOUT, 9012) NPIN 

°012 FORMAT ( 48H I MPROPFR NUMBFR OF NONZERO PRIOR PROBABILITIES 16) 

GO TO 2035 

8031 WR ITEI IUNOUT, 9031) UCOFF 

°033 FORMAT (20H UCOEF OUT OF RANGE G14.6) 

GO TO 2035 

8035 WR ITEI IUNOUT, 9035) TUTILF 

9035 FORMAT ( 31H UTILITY FUNCTION CHOICE ILLEGAL 16) 

GO TO 2075 

8038 WR ITF( IUNOUT, 9018) P 


575 


576 


577 


576 


575 


5 8 C 


501 

814 

582 


503 

816 

5 84 

617 

585 

819 

586 


587 

821 

580 


5 8 9 


590 


591 


592 


593 


5 94 


595 


596 


597 


598 


599 


6CC 


601 


602 


603 


6C4 


605 


606 


607 


608 


609 


6 1 C 


611 


612 


613 


614 


615 


616 


617 


618 


619 


620 


621 


622 

823 

623 


624 


625 

825 

626 


627 


628 

827 

629 


63 C 


631 

829 

632 


633 


634 

831 

635 


636 


637 

833 

638 


639 


64 C 

835 
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9038 

FORMAT ( 3 3H ILLEGAL 

INPUT PRIOR PROBABILITY 

G1 4* 6 ) 

641 



GO TO 2035 



642 


8050 

WR I TF ( IUNOUT, 9050) 



643 

837 

90 50 

FORMAT ( 52H ILLEGAL 

CHARACTER USED TO INPUT 

A GENERATOR FOR DPG ) 

644 



GO TO 20^5 



645 


8040 

WR T TF ( IUNOUT , 9060 ) 



646 

839 

*9060 

FORMAT ( 2 8H GENERATORS NOT INDEPENDENT ) 


647 



GO Tn 2035 



64 8 


8070 

WR T TF ( IUNOUT, 9070) 



649 

841 

9070 

FORMAT ( 45H ILLEGAL 

CHARACTER USED TO INPUT 

BLOCK EFFECT ) 

65C 



GO TO 2035 



651 


8140 

WR I TF ( IUNOUT, 9140) 



652 

843 

9140 

FnRMA T ( 2 OH INVALID 

NOo OF DPGS ) 


653 



GO TO 2035 



654 


8150 

WR ITF ( IUNOUT, 9150) 



655 

845 

9150 

FORMAT ( 26H INVALID 

NO, OF GENFRATORS ) 


656 



GO TO 2035 



657 


81 60 

WR I TF ( IUNOUT, 91 60) 



658 

847 

0160 

FnP MAT ( 27H STOPPING 

i PROB OUT OF RANGE ) 


659 



GO TO 2075 



66C 


8170 

WR ITF ( IUNOUT, 95.70) 



661 

849 

9170 

FORMAT f 29H INVALID 

NOo OF BLOCK EFFECTS ) 


662 



GO TO 2035 



663 


8180 

WR I TF ( IUNOUT, 9580) 



664 

851 

9i 80 

EO°MAT ( 3TH BLOCK EFFECT PROB OUT OF RANGE > 


665 



GO TO 20’F 



666 



END 



667 



SIPFTC L INFR 


SUBROUTINE LINF(T,U) 1 

COMMON /BDATA /POWER SO X ) , I ALPHA ( 9 ) , AL PHA (10 ) , I UN IN, I UN OUT , M ASKK, NEG 2 

COMMON /BT/ NF AC ,NCLA$ S, NN * NDPG , PBA YES , IRUN, PBAYEX , IUT ILFtUTSWCH 3 

COMMON /L TNX/ XNOUT ( 5) tHOLOUT ( 9) t FAC (1 4 ,9 ) 4 

D I M C N $ TON XNUMFR f 9),BLANK(3) ,BL0CKS(4) 5 

OATA( BLOCKS( I) ,1 = 1, A) /6HC0NF0U ,6HNDE 0 W,6HITH BL ,6 HOCKS / 6 

DATA(BLANKU) t T=1 t31 '0777777606060, 0606060777777 t 6H / 7 

nATA(XNUMFR( I) ,T=X,9) / 0605 A OX 77777 7 , 0777777605402, 8 

X 0605403777777, 0777777605404, 0605405777777, 0777777605406 , 9 

X P6054CT777777, 0777777605410, 0605411777777 / 1C 

DATA MA SK /P3 / , ZFR0/6H 0* /, G ME AN/6HG MEAN/ 11 

FOU I V A LENCF (X,IX) , (Y,IY) 12 

C 12 

C *********** ************************************************************ 

IX= T 15 

IMTXoFOoO) WRITE ( IUNOUT, 5005) ZERO,GMEAN,U 16 3 

I F ( TXoFO.O) R F TURN 17 

NF= N F AC 18 

JJ=1 19 

NFl =N F - 1 2C 

no 1 00 J=x ,NF 1 , ? 21 

NF= NF-? 22 

Y=AND( MASK,X) 23 

IY= TY + 1 24 

C 25 

GO TO (20,40, 6C, 50) ,IY 26 

20 XNOUT ( JJ )= BLANM3) 27 

HOtniJT( J )= Bt ANK ( 3 ) 28 

HOLnUT(J + l)= B LANK ( 3 ) 29 

GO TO c 5 30 
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c 31 

40 XNOUT! JJ)= AND ( BLA NK ( 1. ) , XNU PER ( J) ) 32 

HOLOUT(J)= FAC ( l, J ) 33 

HOLOUT ( J 4l )= BLANK ( 3) 34 

GO TO 45 35 

C 36 

40 XNOUT ( J J ) = ANtMBLANK! 2) .XNUMER! J4-1 ) ) 37 

HOLOUTIJ ) = BLANK ( 3 ) 38 

HOL OUT ( J +L ) = F AC ( 1 » J+ 1 ) 35 

GO TO 55 40 

C 41 

BO XNOUT ( JJ ) = AND ( XNUMER ( J) , XNUMER ( J+l J 1 42 

HOL OUT ( J ) = FAC ( 1 , J ) 43 

HOL OUT (J + l )= FAC(1,J+U 44 

C 45 

05 TX=IAR S! ?» X) 46 43 

JJ= JJ+1 47 

100 CONTINUE 48 

I F ( M F ) 500,130,105 49 

105 X= AN DIMASK, X) 5C 

IX= 1X41 51 

GO TO (110,120), IX 52 

HO XNOUT ( JJ ) = BLANK I 3 ) 53 

HOL OU T I N F AC ) = BLANKI3) 54 

GO TO 130 55 

120 XNOUT! JJ)= ANDIBLANKIl ) , XNUMER I NFAC) ) 56 

HOL OUT (NFAC ) = EACH, NFAC) 57 

C 58 

C 59 

130 CONTINUE 6C 

WRITE! IUMOUT, 5000) (XNOUT! I) ,1=1 ,5) .(HOLOUTII) ,1 = 1,9), U 61 6C 

400 RETURN 62 

5000 FORMAT ( !H 5A 6 , 6X , 9A 6, G1 4. 6 ) 63 

5005 epRMATflH A6,30X,A6,43X,G14.6) 64 

5010 FrtRMATOH 40X,A6) 65 

ENTRY ALINE(TA,IALIAS) 66 

DIMENSION IALIAS(l) 67 

00 BOO K =1 , 1 A 6e 

IX=IAL IAS! K) 65 

IF(lx.EQoO) WR I TE ( I UNOUT, 501 0) ZERO 7C 8C 

T F( IXeEO.O) RETURN 71 

NF = NFAC 72 

JJ = 1 73 

N FI = NF-1 74 

00 700 J = 1 ,N F 1 ,2 75 

N F = N F— 2 76 

Y = AND(MASK,X) 77 

I Y= I Y 41 78 

GO TO (620, 640, 660,630), IY 75 

6->0 XNOUT! JJ )=BLANK( 3 ) eC 

GO TO 695 81 

640 XNOUT ( J J ) = AN0(BLANK(1 ) .XNUMER! J)) 82 

GO TO 655 83 

660 XNOUT! JJ) = AN0!BLANK(2) , XNUMER! J + D) 84 

GO Tp 69*= 85 

6 B 0 XNOUT! JJ) = AND(XNUMER(J) ,XNUMER!J4l) ) 86 

69c TX = IAPS(?,X) 87 107 

J J = J J 4 1 88 

700 CONTINUE 85 

IF(NF) 500,730,705 5C 

705 X = AN0 ( M A SK , X ) 91 

ix=rx4i 92 

GO TO (710,720), IX 93 

710 XNOUT! JJ )= BL ANK! 3 ) 94 

Gn TO 73C 55 

770 XNOUT! JJ)= AND(BLANK(1 ) , XNUMER ( NFAC) ) 56 

730 WR ITF! IUNOUT, 850) ( XNOUT! I ) , I =1 , 5) 57 120 
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*00 CONTINUE 98 

P9C FORMAT <40X,5A6) 99 

RFTUPN IOC 

*>00 STOP 101 

END 102 


sibftc ERECT 


SUBROUTINE RFC T(N , X, Y , NFAC ) 

OIMFNSION X(9, 33), Y(33,BB) ,X0UT(9) 

COMMON /PDATA /POWERS PI ) , I ALPHA (9) , ALPHA (10) • I UN IN, IUNOUT , M ASK, NEC 
OATA J 9/ 9 / 

LOGIC A L OUT 
OUT= o FALSFo 
J T IMF S =0 
JCOL= N 
TW = *J 
9 TWW=1 

j NXT=J COL-9 
IF(JNXT) 10,10, BO 
i C JP = JCCL 
OlJT=o TRUE© 

GO TO 50 
JCOl = JNXT 

Jp — J o 

90 GO TO ( 5 2, 94) , IW 
TW= JP- 1 

WR!TF< IUNOUT, 1001) (J,J=1 ,IW) 

WR T T E ( IUNOUT, 1 002 ) 

T W= 7 

GO TO *6 
54 J1 =JTIMFS 
JL - Ji +JP-? 

WR ITF( IUNOUT,! COB) < J , J=J1 ,JL) 

WR ITE ( IUNOUT, ICO?) 

96 OH *OC 1=1, N FAC 
on 70 J =1 , JP 

JJ=JTIVFS+J 
TO XOUTf J ) = XU, JJ) 

WR I tr ( IUNOUT, 1C09) (XOUT(K) ,K=X ,JP) 

100 CONTINUF 

WR I TF ( IUNOUT, 1002) 

OH 200 1*1, N 

nn pc J=X,JP 
JJ= JT1MFS+J 
170 XOUT( J )= YU, JJ) 

GO to (177,1.74), IWW 
177 WR ITF< IUNOUT, 1007) (XOUTfK ) ,K=1 , JP) 

I WW = 7 

WR TTF( IUNOUT, 1008) 

GO TO 200 
174 JI « T-' 

WP ITE ( IUNOUT, 1009) JI , ( XOUT(K) ,K = 1, JP) 

TOO CONTINUE 

I F ( OUT ) PF TURN 
J T IMF S = JT IMF S + JP 
GO TO 9 

1001 F0RMAT(71H!SUMMARY OUTPUT TABLE / ! HK/8 X , ]. 2H** BAYES ***,8(6H *** 

X* 12, 6H ***** ) ) 

1002 FOP MAT ( 1H ) 

1007 F0RMAT(7!H1SUMMARY OUTPUT TABLE /1HK/9(6H ****12, 6H ***♦*)) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 

2C 15 

21 21 

22 

23 

24 

25 

26 25 

27 31 

28 
29 
3C 

31 

32 41 

33 

34 48 

35 

36 

37 

38 

39 

4C 59 

41 

42 65 

43 

44 

45 68 

46 

47 

48 

49 
5C 

51 

52 

53 
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r> r> n r» r> 


' 005 FpRM AT ( 6X * 8( 7X * Al * 6X) ) 54 

1 00 7 FQRMATM 7HKEXPECTED UTILITY /1 8H OVER STOPPING PTS /1HK/ 55 

X 6X,5G14. 5) 56 

100 8 FPPMATO 7HKEXPFCTED UTILITY /20H AT EACH STOPPING PT ) 57 

1008 FORMAT M.H 13, 2X, 9G14o 5> 58 

END 5S 


tTBFTf, PERM DEBUG 

SUBROUTINE PERMUTC IA»N ,LOG , I SEND I 1 

COMMON/BDATA/ POWFRS(ll), IALPHA(O), ALPHAUO), 2 

X I'JNTN, IUNOUT, MASK, NEG 2 

COMMON /BT / NFAC >NCLASS* NN *NDPG » PB AYES * IRUN,P8AYEX , IUT ILF, UTSWCH 4 

COMMON /B2/ DP G ( 1 28 ,32 1 , BLOCK ( 1 28 ,32 ) ,IDPG(32I , NBLOCKI 32) , I P( 9 ) , 5 

X ID( ° ) , N SUB 1(9) ,LPERM(8),I I UO) ,PST0P(32) ,WT<32 > , PSUMX { 32 ) , PSUMC 32 ) 6 

INTEGFR POWFR S,DPG, BLOCK 7 

LOGICAL LOG 8 

EQUIVALENCE ( X , I X ) , ( SI , I S ) , ( S J , J S) 5 

DIMFNSION IA(1),MASK1(15) ,MASK2<15> 1C 

EQUIVALENCE { MASK1 , POWER S ( 2 > ) 11 

DATA! MASK7(I) , 1-1,0 1/07777777774, 07777777771 ,07777777763, 12 

X 07777777747,077777777) 7, 07777777637, 07777777477, 13 

X 07777777)77,07777776377 / 14 

15 

*** * *** *** *4******* ************************************ *** ************* 15 

N T =N 17 

IF(NT-l) eqo ,500,5 18 

5 IFI^NOToLOG) GO to 20 1<3 

DO 10 K=?,NT 20 

IP(K)= 0 21 

T 0 ( K ) - 1 22 

’0 CONTINUF 22 

LOG- oFALSE, 24 

25 

* ** ***** ** *** ** ******** **** ********** ***** **** **** ********** *** ******** 26 

70 K-0 27 

3 C IQ- IP (NT ) + ID(NT) 28 

IP ( NT ) = 10 25 

IF(IQ-NT) 80,40,80 3C 

40 ID( NT ) = -1 31 

45 IF(NT-5)60,60,50 32 

50 NT -NT - 1 33 

GO TO 30 34 

60 TO- 1 35 

LOR- 0 7RIJF 0 36 

GO TO ) 6 C 37 

80 IF( IQ) 1 50,85,150 38 

8 C ID( NT ) =1 35 

K= K + 1. 40 

GO TO 45 41 

42 

****** ********* *********************************** ************* ******** ^2 

1 5 0 TO- TQ+K 44 

I- IA ( IQ ) 45 

I AC TQ ) = IACTQ+lt 46 

IAI TQ + 1 )= I 47 

12= ISFMD+IQ 48 

M) = MA SKI ( 17-1 ) <45 

M2— M A SK 1 1 12) 50 

M3- M A SK 2 ( 12-1) 51 

DO 400 1=1 ,NDPG 52 

ID I = I DPG ( I ) 53 

IB- M BLOCK ( I ) 54 
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rv n n n o 


*************************** 

I F ( mi ) ?40,24C,1 90 

i9o on ?oc m=i, idi 

IX = 0PG(M,T) 

SI= ANC(M1 tlX ) 

SJ= AND ( M2* I X ) 

X = AND (M3, IX ) 

IS= IS*? 

JS=JS / 5 

DPG(M, I)=IX+IS+JS 
poo CONTINUE 

*************************** 
240 TF(IB) 4C0,400,P50 
PRO on 300 M=l, IR 
IX= BLOCK ( M* I ) 

S T = ANC(M’ ,1X1 
SJ = AN0(M2,IX) 

X = AND (HP, IX ) 

IS = IS*? 

JS=JS/P 

BLOCK ( M, I ) = IX+IS+JS 
POO CONTINUE 
400 CONTINUF 
ROC RFTUPN 
END 


IRFTC P.FMAX DERUG 


SUBROUTINE RFMACH 1 

COMMON /BDATA / POWERS(ll), IALPHA(9), ALPHA! 10 1, 2 

X I UN I N , I UN OUT, MASK, NEG 3 

COMMON /B1 / NFAC,NCLASS,NN,NOPG, PBAYES, IRUN.PBAYEX , IUT ILF,UTSWCH A 

COMMON /BP/ 0PG(128,32) ,BLOCKO. 28,32) ,IDPG(32) , NBLOCK( 32 ) , I P( 9 ) , 5 

XTD(°),NSUBI(9) ,LPERM(9),imO) .PST0PC32) , WTO 2 ) , PS UMX { 32 ) , PSUM ( 32 ) 6 

C0MMON/REST/PR0B(512) , TU512), PBLOCK (1. 2 8 ,32 ) , I S AVEP ( 9,33), 7 

X UTIL ('■)?), ULISTU2B), I UL I ST (1. 28 ) , I00CG(4,32) 8 

INTEGER POWERS, T1 , DPG, BLOCK 9 

COMMON /L INX/ XN0UT(5) ,H0L0UT(9) ,FAC(14,9» 10 

(OGICAL UTSWCH ,ONLVI 11 

DIMENSION K2C YCL < 9 ) , KPER M ( 9 ) 12 

DIMENSION M A SKO ( 9 ) ,MASKI. ( 9! ,KKSAVE(9) 13 

fouIVALENCE(MASK0,P0WEP.S(2) ) 14 

DATA (MASKKI ), 1=1,9) / 0T777777776, 07777777775, 07777777773, IS 

X 077 7 777T767, 07777777757, 0T7777 7773 7, 07777777677, 07777777 577, 16 

X 07T77777377 / 17 

EQUIVALENCE ( X , I X ) , ( I S , SI ) , ( JS , S J) 18 

DIMENSION IALIASU28) , SUMALF (9,33) , SUMVAL (33 , 33 ) 19 

20 

*********** *** ************* **** ******** **** ******* ************* ******** 2 1 

22 

R FM AC V- SFCTION \ 22 

C 24 

HNLY1= e FA L SF 0 25 

NDPGl^ NDPG+l 26 

GO TO 4 27 

r NTRY ONCF (* ) 28 

NOPGl = 1 29 

ONf. Y^ = oTRUFo 20 


55 

56 

57 

58 
55 
60 
61 
62 

63 

64 

65 

66 

67 

68 
65 

70 

71 

72 

73 

74 

" 75 

76 

77 

78 

75 
8C 
81 
82 
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.-inn r> o o r> n r*> n o r> nnn d o t n 


A OH 10 IL=1,NFAC 31 

KKS= ISAVFP(LL,1) 32 

KPPRM< KKS)=LL 33 

*> 0 KK S A VE ( KK S 1 = L L 3A 

35 

******* *****************************************************+** ******** 2 6 

37 

on 250C L=1 »N0PG1. 38 

39 

************************************* ********************************** AC 

A 1 

REMACH SECTION 2 42 

43 

IFI0NLY1 ) CO TO 25 44 

NODT= L-1 45 

I F I NO UToEOoO) WRITE(IUNOUT»5030) 46 

I F ( NO It To NE » 0 1 WR I TE (I UNOUT , 5035 ) NOUT , ( I DDCG (K • NOUT ) ,K= 1 , 4 ) 47 

25 wp ITE< IUNOUT, 40901 48 

00 TOO LL =' » NF AC 45 

KKS = I S AVEP ( LL » L) 50 

WRTTF( IIJNOUT, 4095) LL, FAC <1 . LL) , ALPH A< KKS ) 51 

SUMALF (LL,L)= ALPHA (KKS) 52 

TOO CONTINUE 53 

54 

*********************************************************************** 5 5 

56 

R FMAC K SFCTION 3 57 

58 

no 120 LLL=7. tNFAC 55 

K2CYCL (LLL )= KPFRM(LLL) 60 

on 1.T5 K= LLL ,NFAC 61 

IF(KPERM(K)-LLL) 115,112,1.15 62 

112 Kpfrw|!K)= KPFRM(LLL) 63 

GO TO 120 64 

1.1 5 CONTINUE 65 

’20 CONTINUF 66 

67 

**** ******* *********** ***************4************************* ******** 6 8 

65 

REMACF SFCTION 4 70 

71 

N F 1=N F AC- 1 7 2 

00 158 LLL =1 , NF 1 73 

I = NFAC-LLL 74 

J= K’CYCLI I) 75 

KK= J- 1 76 

IFIKK 1 1 58,1.58,130 77 

130 M 0 T = MASKOd 1 7e 

MOJ= MASKO(J) 75 

MU = MASKKI) 80 

Ml J= MS SKI ( J > 81 

82 

****************** 83 

84 

00 157 L4= 1 ,NDPG 85 

14= IOP G( L4 ) 86 

IFI 14 ) 15">, 1.52, 148 87 

14P 00 150 L 5=1 , 14 88 

IX = DPG(L5,L41 85 

S I = AND(M0I,X1 90 

SJ= ANC(MOJ,X) 91 

X= AND(AND(M1 I,X) ,M1J) 52 

IS=TALS(KK,SI ) 93 

JS=IAR$(KK,SJ) 94 

0PG(L5,L4)=IX+IS+JS 95 

i 50 CONTINUE 96 

C 97 


23 

25 

31 


35 


85 

87 


51 


C* ************************** gg 

16? T 4 = N BLOCK( L 4 ) <59 

IF( 14 ) 167,1*57,1 55 IOC 

155 DO 1 56 L 5 = 1 , T 4 101 

I X = BL OCK < L 5, L 4 I 102 

S I = ANOCMOItX) 103 

$J= ANC(MOJ,X) 104 

X- AND (AND CM 1 . I , X ) , M1.J ) 105 

TS=TAL 5fKK,ST ) 1C6 101 

JS=IAR S<KK, SJ ) 107 103 

Bt OCK( L5 f t.4)-IX+IS+JS 10 8 

i ^6 CONTTNLF 105 

i *-> CONTINUE 110 

158 CONTINUF 111 

r 112 

r *********************************************************************** 113 

r 114 

C REMACI- SECTION 5 115 

C 116 

DO iOOC T=1,NDPG 117 

WR ITFt IUNOUT, 5070) I 118 113 

WRTTF( IUNOUT, 5022) ( I DOCG (K,I),K-l»4J 115 114 

I RUN = T RIJN + 1. 120 

PMJ M X ( n=OoO 121 

NB=NBL OCK ( I ) 122 

TOI = IDPG ( I ) 123 

1011=101+1 124 

C 125 

r *-*-*_*_ *-*_*-.*-*_*_*-*_*_*-*-*-*-*_*_*_*_*_*_*_*_ *_*_*_*_*_*_* _*_*_* 126 

r 127 

C 128 

TE(NB) 7010, 728,1 66 129 

166 on ??6 K=1 , NB 130 

KT= Bl. OCK(K,T) 131 

KI1 = K T4l 132 

T 1 ( K 1 1 1= IP UN 133 

1A=0 134 

I E ( IDT ) 701.0, 167,1 71 135 

C 136 

C 137 

F 138 

167 loiJT= K 1 1 13S 

U= MTIMKIl) * ( 1. o 0— P 8 L OC K ( K , I ) ) 140 

OH TO 220 141 

C 142 

C j 4 3 

C 144 

171 IF(PBLOCK(K, I )-1 o 0)176,1 73,7010 145 

173 OH 174 IK= 1 , TOI 146 

JJ= IEXnR(KI,DPG( IK, I ) ) + 1 147 148 

T1 ( JJ )= IRUN 148 

174 CONTINUE 145 

TntjT= K 1 1 150 

U =0 0 151 

GO to 220 152 

C 153 

r 154 

r 155 

176 PS = P R 0 B ( K 1 1 ) 156 

PR= PS 157 

IUL 1ST (1 ) = KT1 158 

UL T ST ( 1 ) = UTILfKIU 155 

TSTAR=C 16C 

14 = 0 161 

I E ( P $ ) 1.78,178,175 162 
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-> n n non non non 


^8 T ST AR = 1 

ultst= -utiliktu 

I At. !AS(1)= KT1 
PR=t o 0 

1 79 IM4X= KIT 


IF(UTSfcCH) GO TO ?00 
nn 187 KK=1, IDT 

JJ=IFXOR(KT,DPG(KK,I) ) +1 
Ti ( JJ ) = TRUN 

IF(PPOB(JJ)1 7010,184,18? 
1.8? PR = P R *PROB( J J ) 

I F ( PR 0 8 ( J J )-P S 1 186,183,187 
183 T A= I A + 1 

TAL I A S f T A ) J 
GO TO 18 7 
T 84 ISTAR= ISTAR+1 

IAL TASUSTAR )=JJ 

186 PS= PFOB(JJ) 

TMAX= JJ 
TA=C 

187 CON TIN IE 


IF( ISTAR-1 ) 190,192,1 89 
i 89 U=0 0 

TA = TSTAP 
imjr= TM AX 
GO TO 220 

i P C PR= PR /P RD B ( IMAX) 

l.o? u=PR*L7IL< IMAX)*(1*0-PBL0CK(K, I ) I 

iout=tmax 

GO TO ??0 


’00 00 206 KK =2, IOIT 

J J=!FXOR C K T ,0PG(KK-1, I ) ) + l 
T1 i JJ>= IRON 
IIJLTST(KK) = JJ 
IF(PROB(JJ)) 7010,204,202 
’0’ PR= PR * PROB < JJ) 

UL T ST ( KK ) = UTIL( JJ) 

GO TO ?C6 
704 I ST AR = ISTAR + 1 

I AL TA S ( ISTAR)=JJ 
TM A X = JJ 

ML IST( KK )= -UTIL( JJ) 

’06 CONTINUE 


I F { ISTAR-1) 708,214,1 89 
208 00 210 KK = 1, T0T1 
JJ= I UL I ST (KK ) 

1JL IST(KK)= ULI ST(KK)*PR/PROB ( JJ) 
2"* 0 CONTINUE 
1JMAX=0 3 0 
IOUT= I UL I ST( 1 ) 

DO 21 3 KK = 1, 1011 

TF(HLIST(KK) - UMAX) 213,212,211 
?H t)MAX= ULIST(KK) 

T A = 0 

I O UT= I UL I S T ( K K ) 

GO TO 21 3 


163 

164 

165 

166 
167 
i6e 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 
189 
19C 

191 

192 

193 

194 

195 

196 

197 

198 

199 
20C 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 
212 

214 

215 

216 
217 
216 
219 
22C 
221 
222 

223 

224 

225 

226 

227 

228 
229 
23C 


172 


2C7 
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r> ^ n non r> o o no r> 


in 


231 

232 

233 
23 ^ 

235 

236 

237 

238 


714 00 2’ 6 KK =1 , 1011 235 

KS = KK 240 

TST = AND(ULIST(KK),NFG) 241 

TST= OR(TST,MASK» 242 

IF! TST ) 218,218,21 6 243 

2’ 6 CONTINUE 244 

”8 iniJT=IlLTST(KS) 245 

II=-IJLI?T(KS!*PR*(1,0-PBL0CK!K,I) ) 246 

247 

248 

245 

770 psumxi n=PSUMxm+u 25c 

toot = inuT-i 251 

CALL L INF ( IOUT.tj) 252 272 

IF! IA ) 726,726,22’ 253 

7’2 00 223 1 1 A= 1 , IA 254 

22’ IAL TAS(TIA1 = TAL!ASUTA)-1 255 

CALL A L INE (IA, I ALIAS) 256 282 

2’6 CONTINUF 257 

258 

-*-*-*-*-*-*~ *-*-*-*-*-*-*_*-*-*-*-*-*-* -*-*-*-*-*-*-*-*-*-*-* -*-*-* 255 

260 

228 TF(ini) 7010,225,245 261 

***************** 262 

225 00 240 K=1,NN 267 

T F( T1 < K I- IRON ) 230,740,730 264 

730 U=U T IL (K. I 265 

’35 T0U7= K— 1 266 

DSUMX(n=PSUMxm+U 267 

CAl L L IN P ( IOUT , LI I 268 58 

’40 CONTINUE 265 

GO tq 1000 270 

271 

272 

273 

’4B no TOC K=1,NN 274 

I F( T7 ( K I-IRUN ) 246,700,’46 275 

746 PS= PROB(K) 276 

KM7 = K-7 277 

PR = P S 27 8 

TST AR = C 275 

T A = 0 280 

IOL 1ST!’ 1 = K 281 

ULISTI 1)= UTIL (K) 282 

TP(PS) 248,243,750 287 

’48 T8TAR=7 284 

I A! TA S ( I )= K 285 

ULTST! 71= -UTIL(K) 286 

dr=i,0 287 

’50 TMAX=K 288 

IF(IJTS7>CH) GO TO 430 285 

C 250 

00 4T 0 KK =1,101 251 

J J=IEXOR(KM’ ,OPG(KK,I ) )+l 232 325 

T!(JJ)=TRtJN 253 

TP(PRDB(Jjn 380,330,380 284 

’80 PR = PR TPROB! J J ) 255 

IFIPROB! JJI-PS) 400,385,410 256 


712 I A= I A + 1 

IALIAS(IA)=IULIST(KK) 

2 ' 3 CONTINUE 

U=HMAX*( i ,0-PBL0CK(K, I ) ) 
GO TO 220 
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n n n n n n n on on n 


5 

I A= TA +1 

297 


IAl. I A S f IA )=JJ 

298 


GO Tn 410 

299 

^90 

ISTAR = TSTAR + 1. 

200 


IAl IASC I STAR )=JJ 

301 

400 

PS=PRDR( JJ) 

302 


TMAX=J J 

302 


!A = 0 

304 

4’ 0 

CONTINUE 

305 


GO Tn 50*5 

306 



307 






309 

430 

on 460 KK=2, ID 11 

310 


JJ=IFX0R(KM1 ,DPG(KK-1, I) )+l 

311 


TM JJ >= TRW 

312 


ItlL I ST (KK ) = JJ 

313 


TF(PRPBtJJ)) 450,450,440 

314 

440 

RR= PR *PRnB( J J ) 

315 


IJL I ST ( KK )- UT I L ( J J » 

316 


GO Tn 460 

317 

* 50 

T STAR = T STAR + 1 

318 


TALTASI I STAR )=JJ 

319 


TM AX = JJ 

320 


ULTSTf KK) = -UTItC JJ) 

321 

460 

continue 

322 


GO Tn 60 0 

323 



324 






326 

*50 * 

TF( ISTAR-1 ) 520,430,51 0 

327 

S 1 0 

U = OoO 

328 


TA=ISTAR-l 

329 


TOUT = I M A X 

330 


GO Tn 690 

331 

5?o 

PR =PR / PRHB ( IMAX) 

332 

C7 0 

U=PR*UTTL( IMAX) 

333 


IOUT= I MAX 

334 


GO Tn 690 

235 



336 






338 

600 

r FC IST AP-I ) 620, 66 C , 51 0 

339 

6?0 

DO 6^0 K K = 1 , IDI1 

340 


JJ= T UL I ST (KK ) 

241 


UL I ST ( KK ) = ULIST(KK) * PR/PROB(JJ) 

342 

6^0 

CONTINUE 

343 


IJM AX = 0«,0 

244 


in IJT = I LL T ST ( 1 ) 

345 


on 650 KK = 1 , mil. 

246 


IFUJl I ST ( KK ) - UMAX ) 65 0,640,635 

347 

69* 

UMAX = ULIST(KK) 

34e 


T A = 0 

349 


inUTssIU. IST(KK ) 

35C 


GO Tn 650 

351 

640 

I A= I A + 1 

252 


IAL IASC IA )=IUL IST(KK) 

353 

6*0 

CONTINUE 

354 


U=!JMAX 

355 


GO Tn 690 

356 



357 






359 

640 

DO 670 KK=1,inil 

360 


KS= KK 

361 


TST = AN H( UL 1ST ( KK ) , NFG ) 

262 


TST=OR (TST, MASK) 

363 


TP(TST) 680,630,670 

364 



351 
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6 T 0 CONTINUE 365 

*RC TOUT=IULIST(K$) 366 

U=-UL I ST(KS)*PR 367 

r, 368 

C 365 

C 37C 

600 PSUMXI I)=PSIJMXm+U 371 

T0UT=ICUT-1 372 

CALL L INF ( TOUT , U ) 373 426 

I F( TA ) 700,700,605 374 

605 00 606 I TA=1 » IA 375 

606 I AL IAS ( 1 1 A 1=1 ALI AS ( II A )-l 376 

CALI AL INF (IA, IALIAS1 377 436 

TOO CONTINUE 378 

tCOO CONTTNUF 379 

C 38C 

C ********** ************************************************************* 381 

C 382 

C REMACE SECTION 6 383 

r 384 

PR AYEX = OoO 3e5 

00 I 0 7 C 1 = 1 , NOPG 386 

PSUMXI 11 = PSUHXm*WT(II 387 

PRA YEX = PB AYEX + P SUM X (I) *PSTOP ( I) 388 

SUMVAL ( I +1 »L ) = PSUMXU) 389 

I.OPO CONTINUE 300 

SUMVALH,LI = PBAYEX 391 

WR ITF( IUNOUT, 5000) PBAYEX 392 453 

WR I TE ( I UNO U T , 5005 ) ( I , PSUMX (I) , I =1 , NDPG) 303 

C 394 

C *************************************************************** ******** 395 

C 396 

C RFMACE SFCTION 7 A 397 

c 398 454 

TF( L-NDPG1 ) 105C, 2700,7700 399 

1 050 00 1 07 C LL=! , NFAC 400 

J = T5AVEP (LL,L + 1 ) 401 

KPERM ( J )=LL 402 

10T0 CONTINUE 403 

c 404 

r* ********** ************************************************************ 60 5 

C 406 

C REMACE SFCTION 7B 407 

00 ?00C LL=). , NFAC 40e 

no i 0 B C L4=’ , N P AC 409 

IF(KKSAVP(L4 )-LL) ’ OBO C75 ,1 080 410 

i 075 K7CYCL (LL)= KPERMI L4) 411 

GO TO 2000 412 

1 0*0 CONTINUE 413 

?000 CONTINUE 414 

C 415 

C *** *********************************************** ****** ******* ******** 416 

r 417 

C REMACE SECTION 7C 418 

00 P01 C LL=l , NFAC 419 

KKSAVE(LL)=KPERM(LL) 420 

KPFRM( LL )= K2C YCL(LL) 421 

20’ 0 CONTINUF 422 

C 423 

f ************************************************** ************* ******** 424 

C 425 

I'OC CONTTNUF 426 

?700 IF(ONLYl) PFTURN1 427 

CALL RECT(NOPG!,SUMALF,SUMVAL,NFAC) 428 496 

P F T UR N 429 

r 430 

r *********** ************************************************************ 431 
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n n o 


■5000 FORMAT ( 5 2HL FOR THF ABOVE PERMUTATION THE EXPECTED UTILITY IS 432 

X 014.5) 43 3 

■5005 FORMAT! 53HK THE EXPECTED UTILITIES AT THE STOPPING POINTS ARE.. / 434 

X ( 2ft H DEFINING PARAMETER GROUP I 3 ,3X,G14.5 ) ) 435 

5010 FORMAT (16H ERROR IN OUTPUT) 436 

50?C FORMAT ( 30HKDEF INING PARAMETER GROUP NO. 13) 431 

502? FORMAT (1.H 4A6/60! 2H *)) 438 

5025 FOPMAT(65(2H ->) 439 

5030 FORMAT ( 37H1 THIS MATCHING IS THE BAYES MATCHING ) 440 

5035 FORMAT ( 5 2H1 THIS MATCHING MAXIMIZES THE EXPECTEO VALUE AT THE 12, 441 

X lf,H STOPPING POINT 4A6 ) 442 

4095 FOR MAT ( 1 H 1 4, 1 X , A 6 ,19X , A1 ) 443 

40«>0 FORMAT (38HKVAR IABL E SHOULD BE CALLED ) 444 

* ** *************** ******** ******************************** ************* 44 5 

446 

FRROR MESSAGF 441 

701 0 WRITF! IUNOUT, R010) 448 458 

3010 FORMAT ( 52H PROGRAM ERROR — PREVIOUSLY GOOD DATA HAS BECOME BAD ) 445 

STOP 450 

END 451 


57 



APPENDIX C 


PROGRAM SYMBOLS 


This appendix presents a listing of the major program variables used in NAMER. 
Dimensioned variables have their dimensions specified. 


ALPHA(IO) 

BLOCK(128, 32) 

DPG(128, 32) 
DUMP1, DUMP2 
FAC(14,9) 
HOLOUT(9) 

IALIAS(128) 


IALPHA(9) 

ID(9) 

IDDCG(4, 32) 

IDENT(14) 

IDPG(32) 

n(io) 

IMAX, IOUT 
IP(9) 

IRUN 

ISAVEP(9, 33) 
ISTAR 

ITYPRN 


First nine letters of the alphabet (excluding I) and a blank; used 
for output of optimal matchings. 

Standard- order subscripts of representative members of alias 
sets confounded with blocks. 

Standard- order subscripts of parameters in d.p.g. 's. 

See section 5 of SPECIAL LEWIS RESEARCH CENTER ROUTINES. 

Hollerith identification of factors. 

Temporary storage for Hollerith output; used in LINER and initial- 
ized in NAMER. 

Standard- order number subscripts of parameters in an alias set, 
besides IOUT, which also maximize expected utility of eq. (6) 
or eq. (7). 

The integers 1 to 9; these are permuted by PERMUT and used to 
indicate optimal matchings for output. 

Indicator vector used in PERMUT. 

Hollerith identification of stopping points. 

Hollerith identification of current problem. 

Number of elements in each d.p.g. (not including identity). 

Subscripts of beginning locations of classes of factors. 

Subscript of parameter which maximizes eq. (6) or eq. (7). 

Indicator vector used in PERMUT. 

Running counter used as indicator of parameters evaluated. 

Optimal matchings. 

Counter for number of parameters in an alias set with prior pro- 
bability 1. 0 of being nonzero. 

Type of run; see section 3 of INPUT DESCRIPTION. 
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IULIST(128) 

IUNIN 

IUNOUT 

IUTILF 

KKSAVE(9) 

KPERM(9) 

K2CYCL(9) 

LD1, LD2 

LPERM(9) 

NBLOCK(32) 

NCLASS 

NDPG 

NFAC 

NN 

NSUBI(9) 

ONLY1 

P 

PBAYES 
PBAYEX 
PBLOCK(128, 32) 
PR 

PROB(512) 

PSTOP(32) 

PSUM(32) 

PSUMX(32) 
SUMALF(9, 33) 


Vector of subscripts of parameters in an alias set. 

Variable input unit designation. 

Variable output unit designation. 

Choice of function; see section 7 of INPUT DESCRIPTION. 

Saves permutation vector required to achieve optimal matching. 
Permutation vector. 

Transposition vector for permutation in KPERM. 

Delimiters for DUMP1 and DUMP2 vectors; see section (5) of 
SPECIAL LEWIS RESEARCH CENTER ROUTINES. 

Logical variables controlling calls to PERMUT. 

Number of alias sets confounded with blocks for each d. p. g. 

Number of classes of factors; see section 3 of INPUT DESCRIP- 
TION. 

Number of d.p.g. 's. 

Number of factors. 

2 nfac 

Number of factors per class. 

Logical variable set to .TRUE, if only a single specified match- 
ing is to be evaluated. 

Temporary storage of input prior probability. 

Overall expected utility of best matching evaluated so far. 

Overall expected utility of current matching. 

Prior probabilities of block effects being nonzero. 

Temporary storage used for calculation of ~j |~(1 - Pj). 

Vector of (1 - p^). 

Probabilities of stopping exactly at each stopping point. 

Expected utility at each stopping point of the best matching for that 
stopping point found so far. 

Expected utility at each stopping point of the current matching. 
Saves optimal orderings for output of summary table . 
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SUMVAL(33, 33) 

TMAX 

TMAXX 

Tl(512) 

UCOEF 

ULIST(128) 

UMAX 

UT 

UTIL(512) 

UTSWCH 

WT(32) 

XNOUT(5) 


Saves expected utilities of various optimal orderings for output 
of summary table. 

Maximum total running time permitted. 

Maximum running time for current case. 

Indicator array used in finding all distinct alias sets. 

Constant used to define utility function 5; see section 7 of INPUT 
DESCRIPTION. 

Vector of utilities corresponding to choices of parameters indi- 
cated in vector IULIST. 

Maximum of ULIST. 

Temporary storage used in input of utilities. 

Vector of u^. 

Logical variable used to indicate whether utility function 1 or 2 
or utility function 3, 4, or 5 is being used. 

Weighting values for the stopping points; see section 7 of INPUT 
DESCRIPTION. 

Temporary storage used in LINER for output of numerical identi- 
fication of parameters chosen to be estimated. 
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TABLE n. - ESTIMATE OF TIME REQUIRED BY PROGRAM - BASED ON 


RUNNING OF SAMPLE PROBLEMS 



2 Total time to evaluate all 0. 02 0. 23 3.90 60. 32 


permutations, min 

Time required to print 0. 04 0. 07 0. 19 0. 44 

out results, min 

Number of d. p. g. ’s 4 4 5 5 

Time to evaluate all 0.005 0.058 0.780 12.06 a 151.0 

permutations divided 
by number of d. p. g. T s, 
min 

3 Total time to evaluate all 0. 03 0. 27 4. 64 75. 67 

permutations, min 

Time required to print 0. 04 0. 08 0. 22 0. 47 

out results, min 

Number of d. p. g. ’s 4 4 5 5 

Time to evaluate all per- 0. 008 0. 068 0. 928 15. 13 a 224. 0 

mutations divided by 
number of d. p. g. T s 

a Estimated from fig. 2. 
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8 

7 

6A 

6 

5 

4 

3A 

3 

2 

1 


►/ Subscript Probability Utility 

/ (911) ( FlOO ) (F10.0) 

^ / Number nonzero utility Constant 

( probabilities (16) function (16) (F10. 9) 

| Number in each class (916) 

Number of classes (16) 

Factor identifications (13A6, A2) 

■y Number of factors (16) 

— ( Single specified matching (9A1) * 

Type of run (16) 

Maximum running time (F6.0) 
f Identification (13A6, A2) 


Figure 1. - Pictorial representation of data card arrangement (Asterisk denotes cards which are optional. 
Presence depends on input information contained on earlier cards. ) 
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Number of factors, 
NFAC 

Figure 2. - Logarithm 
of average time re- 
quired per defining 
parameter group as 
function of number 
of factors and choice 
of utility function. 
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